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Abstract 

✓ 

This  report  investigates  the  performance  of  an  active  sonar,  which  is 
limited  by  reverberation,  although  noise  and  limitations  due  to  both  noise 
and  reverberation  are  also  discussed.  If  the  transmitted  pulse  is 
narrowband,  the  difference  between  the  frequencies  scattered  by  the 
target  and  the  reverberation  due  to  their  (hopefully)  different  motion 
relative  to  the  sonar  is  used  to  reduce  the  effect  of  reverberation.  If  the 
transmitted  pulse  is  broadband,  the  ability  of  the  processor  to  decrease 
the  response  to  reverberation  emanating  from  scatterers  whose  radial 
distance  to  the  sonar  is  different  from  the  radial  distance  to  the  target  is 
used  to  reduce  the  effect  of  the  reverberation.  In  either  case,  the 
improvement  in  performance  against  reverberation  due  to  the  processing 
is  calculated  using  various  assumptions  on  the  nature  of  the  reverberation 
and  the  processing.  The  effect  on  the  results  of  varying  many  of  these 
assumptions  is  also  investigated.  The  results  may  be  used  to  investigate 
under  what  circumstances  a  sonar  is  limited  by  noise  only,  by 
reverberation  only  or  by  both.  In  any  case,  the  performance  of  the  sonar 
given  the  average  signal  power,  average  reverberation  power  and  noise 
spectrum  level  and  other  quantities  related  to  the  assumptions  made  may 
be  estimated. - 
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1.  Introduction 

An  active  sonar  system  aims  to  detect  the  reflected  sound  from  a 
target  against  a  background  of  interference,  in  estimating  the  performance 
of  the  active  sonar,  it  is  desirable  to  be  able  to  calculate  the  improvement 
in  signal  to  interference  ratio,  which  is  brought  about  by  the  processing  of 
the  output  from  the  beamformer.  This  report  is  concerned  with  this 
calculation  when  the  interference  consists  of  white  noise  (ie  with  a  flat 
spectrum)  or  reverberation  (ie  sound  scattered  from  all  scatterers  but  the 
target)  when  using  either  of  two  types  of  processing  which  are  described 
in  the  next  two  paragraphs.  Some  of  the  results  in  this  report  are  available 
elsewhere,  but  not,  to  the  author's  knowledge,  in  a  form  which  provides 
sonar  modellers  with  both  the  necessary  basic  results  and  the  mechanism 
to  assess  the  effect  of  changes  to  some  important  parameters. 

One  processing  option,  called  doppler  processing  here,  exploits  the 
fact  that  the  reverberation  return  and  the  target  return  may  be  at 
different  frequencies  since  differences  in  the  motion  of  the  scatterers  and 
the  motion  of  the  target  relative  to  the  sonar  may  induce  different  doppler 
shifts  in  the  returns.  The  output  of  the  beamformer  for  a  time  interval  of 
length  T  (the  analysis  time)  is  divided  into  frequency  bins  by  convolving 
the  output  with  a  (possibly  shaded)  CW  pulse  of  the  desired  frequency.  If 
the  doppler  shift  induced  by  the  target  is  sufficiently  different  from  any 
shift  induced  by  most  of  the  scatterers,  there  will  be  little  or  no 
reverberation  in  the  bin(s)  containing  the  target  signal  and  detection  of 
the  target  will  be  limited  by  noise  alone.  This  report  will  consider  how  the 
performance  of  the  sonar  depends  on  the  difference  in  the  doppler  shifts. 
The  crossover  from  reverberation  limited  to  noise  limited  performance  of 
the  sonar  may  then  be  deduced.  Several  assumptions  regarding  the  details 
of  this  processing  option  and  several  different  models  of  reverberation 
will  be  considered. 

The  above  approach  to  reducing  the  effect  of  reverberation  relies  on 
using  a  narrowband  (eg  CW)  transmitted  pulse.  An  alternative  approach  to 
reduce  the  effects  of  reverberation  is  to  use  a  wide  band  transmitted 
pulse,  such  as  a  chirp  pulse  where  the  frequency  increases  at  a  constant 
rate  for  the  duration  of  the  pulse.  By  convolving  the  signal  from  the 
beamformer  with  a  pulse  with  the  same  frequency  behaviour  as  the 
transmitted  pulse,  the  contribution  from  all  scatterers  outside  an  annulus 
at  a  particular  range  and  of  width  about  c/2B,  where  c  is  the  speed  of 
sound  and  B  the  bandwidth  of  the  pulse,  is  reduced  significantly.  When  this 
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particular  range  coincides  with  the  distance  to  the  target,  the  signal  to 
reverberation  power  ratio  is  increased  due  to  the  reduced  effective 
scattering  area. 

In  section  2,  standard  expressions  for  probability  of  detection  and 
false  alarm  are  given.  It  is  noted  that  under  certain  assumptions,  for  a 
given  probability  of  false  alarm,  the  probability  of  detection  depends  only 
on  a  single  quantity  (which  is  half  the  output  signal  to  interference  ratio 
and  is  denoted  P  and  defined  by  equation  [2.9])  and  the  relationship  is 
shown  in  figure  1 . 

The  quantity  P  is  calculated  in  section  3  when  the  interference 
consists  of  white  gaussian  noise  only  in  the  two  cases  of  a  CW  pulse  and 
an  FM  chirp  pulse  and  well  known  expressions  for  the  detection  thresholds 
are  derived  in  these  cases. 

Section  4  begins  with  a  discussion  of  reverberation.  A  base  case  is 
established  from  which  a  comparison  of  processing  options  may  be 
established.  The  ‘gain"  when  a  processing  option  is  used  is  defined  as  the 
ratio  of  the  value  of  P  when  the  option  is  used  to  the  value  of  P  in  the  base 
case  (the  latter  value  of  P  being  simply  the  average  signal  to 
reverberation  power  ratio).  Fordoppler  processing,  the  "gain"  is  critically 
dependent  on  the  assumptions  made  about  the  nature  of  reverberation  and 
various  forms  of  the  reverberation  amplitude  covariance  function  are 
considered.  Although  ideally  the  reverberation  will  be  concentrated  at  the 
pulse  frequency,  it  may  be  distributed  in  frequency  due  to  motion  of  the 
scatterers  or  uncompensated  motion  of  the  sonar.  Graphs  are  presented  in 
section  4.1  of  the  relationship  between  the  "gain"  and  the  difference 
between  the  doppler  shift  induced  by  the  target  and  the  central  doppler 
shift  induced  by  the  other  scatterers  (suitably  non-dimensionalised). 
Different  assumptions  of  reverberation  frequency  spread  and  different 
combinstions  of  functions  to  shape  the  transmitted  pulse  and  shade  the 
comparison  CW  pulsewith  with  which  the  return  is  convolved  are 
considered.  For  the  FM  chirp  pulse,  the  "gain*  is  not  dependent  on  the 
frequency  spread  of  the  reverberation  and  a  graph  of  "gain"  (figure  22)  as  a 
function  of  bandwidth  is  presented  for  comparison  in  section  4.2. 

In  section  5,  the  nature  of  the  graphs  with  large  differences  in 
doppler  shift  (doppler  processing)  or  large  bandwidth  (FM  chirp  pulses)  are 
explained  by  calculating  asymptotic  series  for  the  gain  in  each  case.  These 
series  provide  a  means  of  estimating  the  "gain*  (relatively  easily)  in 
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situations  not  covered  by  the  graphs. 

Section  6  provides  an  example  illustrating  the  results  of  this  report 
and  section  7  provides  a  summary  of  the  results  presented  in  this  report. 

In  annex  A,  expressions  for  probability  of  detection  and  false  alarm 
are  derived  under  more  general  conditions  than  used  in  the  main  body  of 
the  report  and  these  expressions  are  shown  to  reduce  to  the  standard 
equations  under  the  more  restrictive  assumptions.  These  more  restrictive 
assumptions  are  listed  in  annex  B,  where  the  effect  of  varying  each  of  the 
assumptions  is  also  noted.  Assumptions  made  in  deriving  the  results  of 
sections  4  and  5  are  listed  in  annex  C  and  effects  of  varying  each  of  these 
assumptions  are  described. 
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2.  The  Probability  of  Detection  and  False  Alarm 
A  pulse  of  the  form 

z(t)  -  (2J)0-5  W(t)  cos  (  M(t)  +  2nFt )  [2.1] 

is  transmitted  for  a  time  K.  The  function  W(t)  defines  the  pulse  shape  and 
is  zero  unless  0£t£K  and  is  normalised  so  that 

rK 

K‘1j0  W2(t)dt  -1.  [2.2] 

F  is  the  centre  frequency  of  the  transmitted  pulse  and  M(t)  defines  the 
frequency  modulation  which  is  small  compared  with  F  i.e.  M'(t)/F  «  1 .  The 
quantity  J  is  the  transmitted  pulse  average  power  level  (if  terms  of  order 
(2itFK)'1  may  be  neglected).  Subsequently,  the  sound  received  by  an  array 
of  hydrophones  consists  of  background  noise,  reverberation  and  the  sound 
scatterered  from  the  target  (if  present).  On  the  basis  of  this  sound 
received  by  an  array  of  hydrophones  a  decision  is  made  as  to  whether  the 
target  is  present  or  not.  The  sound  at  the  individual  hydrophones  is  first 
combined  into  beams,  each  of  which  contains  sound  coming  mainly  from  a 
particular  direction.  Each  beam  output  is  monitored  for  a  finite  time  and  is 
often  monitored  for  many  possibly  overlapping  short  time  intervals. 

Suppose  that  one  of  these  intervals  begins  a  time  -cSK  after  the  start  of  the 
pulse  transmission  and  redefine  the  origin  of  time  t  to  be  the  start  of  the 
analysis  period  which  is  of  length  T.  Let  y(t)  be  the  output  of  a  beam  from 
the  receiver  hydrophone  array  for  OstsT  and  on  the  basis  of  y(t)  it  must  be 
decided  whether  there  is  a  target  present  in  the  beam.  It  is  assumed  that 
this  is  done  by  evaluating  the  quantity  r,  where 


r .  (B2+C2)0-5 

[2.3] 

B  -  }  y(t)  H(t)  cos  (N(t)  +  2nft)  dt 

0 

[2.4] 

C-  f  y(t)  H(t)  sin  (N(t)  +  2jrft)  dt. 

fi 

[2.5] 

The  optimal  functions  H(t)  and  N(t)  to  be  used  depend  on  the  properties  of 
the  pulse  and  the  expected  interference.  Usually  N(t)»M(t)  as  the 
contribution  to  the  value  of  r  from  the  target  echo  will  then  be  maximised. 
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The  function  H(t)  is  normalised  so  that 

T'1  J0H2(t)dt  =  1.  [2.6] 


Values  of  r  may  be  calculated  at  more  than  one  value  of  f  near  to  the 
transmitter  frequency  F.  A  detection  will  be  declared  if  r  exceeds  some 
threshold  R. 

Suppose  that  a  coherent  target  echo  of  the  form 

s(t)  -  (2P)0-5 W(t-D)  cos  (  M(t-D)  +  2jtf,t  +0,)  OsD<t£EST 

[2.7] 


is  present  during  part  or  all  of  the  analysis  period,  where  P  is  the  average 
signal  power  level  in  the  beam  (neglecting  terms  of  order  (27rf1K)‘1),  D  is 

the  time  of  the  start  of  the  target  echo,  E  -  min  (T,D+K)  is  the  time  of  the 
end  of  the  target  echo  within  the  analysis  time,  f,  is  the  centre  frequency 

of  the  signal  (approximately  the  centre  transmission  frequency  F)  and  01  is 


the  signal  phase.  It  is  shown  in  Annex  A  that,  provided  certain  assumptions 
are  valid,  the  probability  of  detecting  the  signal  in  a  single  look  is 


P<j" 


e-(u+p)  i0(2V(uP))  du 


-|n<Pfa) 


[2.8] 


where  p!a  is  the  probability  of  a  false  alarm  due  to  the  interference  level 
exceeding  the  threshold  R,  lQ  is  the  modified  Bessel  function  of  zero  order 
and  P,  which  is  half  the  signal  to  interference  ratio  at  the  output  of  the 
processor,  is  give  by 


P-  P  Q2I  (4  a2) 


[2.9] 


where 

Q,2  -  ( H(t)  W(t-D)  cos  (  M(t-D)  -  N(t) )  dt  )2  + 

[E 

( JD  H(t)  W(t-D)  sin  (  M(t-D)  -  N(t) )  dt  )2  [2.10] 
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a  -  J  JT  E[  x(t)x(t-v)  ]  H(t)H(t-v)  cos  (N(t)-N(t-v)+2jrfv)  dt  dv  [2.1 1] 

0  v 

with  x(t)  being  the  component  of  y(t)  due  to  interference  alone  and  E[  ] 
denoting  expected  value.  The  probability  of  detection  as  a  function  of  P  (in 
decibels)  for  various  values  of  probability  of  false  alarm,  based  on 
equation  [2.8],  is  presented  in  figure  1 .  The  assumptions  made  in  deriving 
equation  [2.8]  and  the  effect  of  relaxing  some  of  these  assumptions  are 
described  in  Annex  B. 
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3.  BgafljQQ  in  the  Presence  of  White  Noise 

If  the  interference  consists  entirely  of  whits  notes  with  spectrum 
level  N0,  than 

E[x(t)x(u)]-  N08(t-u)  [3.1] 

where  50  is  the  Dirac  delta  function. 

For  a  CW  pulse  where  M(t)  -  0, 

Q,2  “  (^  H(t)  W(t-D)  COS  N(t)  dt  )2  ♦  (  ^  H(t)  W(t-D)  sin  N(t)  dt  )2 

[3.2] 

fT 

o2-  '<2  N0  o  H*(t)  dt  [3.3] 

so  that 

P  {( !0  H(t)  W(t-D)  cos  N(t)  dt)2  ♦  ( £  H(t)  W(t-D)  sin  N(t)  dt)2} 

P - jj - 

2N0i0H2tDdt  P.4] 

which  has  a  maximum  value  of  PT/(2N„)  when  H(t)  equals  W(t),  N(t)  is  a 

constant  and  the  analysis  period  [0,T]  coincides  with  the  signal  period 
[D.E].  This  is  the  same  as  the  detection  index  (apart  from  a  factor  of  4), 
defined  .for  example  in  reference  [1],  chapter  12.  and  leads  to  the  same 
value  of  detection  threshold.  In  the  case  where  the  signal  and  analysis 
periods  do  not  coincide,  figure  2  shows  the  relative  degradation  in  P  (ie 
decrease  in  P  expressed  in  decibels  from  this  optimal  value)  for  values  of 
D/T  between  0  and  1 ,  assuming  that  E-T  The  different  lines  correspond  to 
different  combinations  of  weighting  functions  W  and  shading  functions  H 
which  are  considered  in  the  next  section  and  are  listed  in  tables  2  and  3 
There  is  a  non-zero  degradation  when  the  shading  and  shaping  functions 
are  different  even  when  D-0.  If  the  analysis  periods  do  not  overlap,  the 
maximum  degradation  is  obtained  when  D/T-0.5  as  for  larger  values  of 
D/T,  the  next  analysis  interval  will  contain  more  of  the  signal  power  and 
will  give  a  degradation  corresponding  to  (1  -D/T)  with  the  weighting  and 
shading  functions  interchanged.  This  maximum  degradation  is  quite  large 
tor  seme  shading  and  shaping  function  combinations  (eg  20  dB  tor  exp 
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shading  and  shaping)  but  can  be  reduced  by  using  overlapping  analysis 
intervals  although  this  increases  the  processing  required  and  may  also 
increase  the  false  alarm  rate  for  a  particular  probability  of  false  alarm 
for  each  analysis  period. 


For  an  FM  pulse, 
fE 

Q,2  *  ( J0  H(t)  W(t-D)  cos  (M(t-D)-N(t))  dt  )2  ♦ 

,E 

( J0  H(t)  W(t-D)  sin  (M(t-D)-N(t))  dt  )2  [3.5] 

J 

a2  -  ’'2No  o^iOdt  [3.6] 


so  that 


P{(^J  H<l)  W(t-D)  OM  (M(t-D)-N(t))  <*)2+(  q  H<t)  W(l-D)  an  (M(l-D)-N(t))  dt)2) 

P - - 

2  N0  J0  H*(t)  dt  [3  7] 

This  is  maximised  when  the  analysis  period  [0,T]  and  the  signal  period  [D.E] 
coincide.  N(t)  differs  from  M(t)  by  a  constant  and  H(t)  equals  W(t).  The 
maximum  value  is  the  same  as  in  the  CW  case. 
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4.  Ptttcfon  in  tfr  Prtttnoi  of 

Let  t  be  a  time  within  the  analysis  penod  [O.T],  Tbs  reverberation  at 
tims  t  comas  from  all  scattarers  whsra  tha  sum  of  ths  travel  timas  from 
ths  transmlttar  to  tha  scattarar  and  from  tha  scattarar  to  tha  racaivar 
lias  within  tha  mgs  (t+t-K.t+t],  where  ths  transmission  startsd  at  tima  -t 
and  tha  puisa  langth  is  K.  Tha  oomponant  of  tha  ravarbaration  from  ail 
scansrsrs  whosa  total  travel  tima  is  C  ♦  t  and  whosa  cantra  frequency  and 
phasa  are  f  and  V  respectively,  may  ba  written 

x(t)  - g (ter.*)  W(t-0 cos (  M(t-C)  ♦  2*ft  +  «’)  t-Ks^st  [4  i ] 

where  tha  phasa  shift  tf  (OSS’S*)  is  due  both  to  scattering  and  wave  travai 
Tha  function  g  which  may  bo  positive  or  negative  includes  aH  factors 
influencing  the  amplitude  such  as  beam  patterns,  propagation  loss, 
scattarar  reflection  coefficients,  interference  between  ravarbaration 
from  different  scattorsrs  and  tha  possibility  of  a  scattarar  changing 
position,  velocity  or  reflection  coefficient  tor  tha  duration  of  tha  puisa. 

Tha  pulse  shape  function  W(t-Q  is  zero  unless  t  s  St-K  Tha 
ravarbaration  from  all  scattorsrs  may  than  ba  expressed  as 

*(t)  -  '0  L  g(t.;\r,m  W(K’)  COS  ( M(t<)  ♦  am  ♦  e1 )  <*;■  dr  <* 

[4.2) 


In  order  to  use  tha  aquations  [A.33]  and  [A.60]  tor  probability  of  detection 
and  falsa  alarm,  it  will  ba  assumed  that  tha  values  of  B  and  C.  constructed 
by  using  x(t)  from  aquation  [4.2]  in  place  of  y(t)  in  aquations  [2.4]  and  [2.5]. 
are  normally  distributed  and  have  zero  mean.  This  assumption  is  valid  if 
there  are  a  sufficiently  large  number  of  independent  scatterers 
contributing  to  tha  ravarbaration.  Tha  expected  value  of  x(t)x(u)  which  is 
required  to  calculate  tha  probability  of  falsa  alarm. 

E[  x(t)x(u)  ]  -  /  j"  M’  M*E[  gac.f.e’)  g(u.;-,r.e-)  jW(t-C)  W(u-;*) 
t-K  u-K  —  —  0  0 

cos  (  M(t-C)-*-2*rt+e’)  cos  (  M(u^Wu^)  de”  d*  df  df  <*'  dC 

[4  3] 


Insufficient  data  are  available  to  estimate  with  confidence  the  covariance 
of  g  which  occurs  in  this  integral.  In  this  report  various  assumptions  are 
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mad*  about  tha  nature  of  this  covarianca  function.  A  spaa  fie  form  of  the 
function  is  considered  in  detail  in  sections  4.1  and  4.2  where  the 
improvements  to  be  made  using  doppier  processing  or  FM  chirp  pulses 
respectively  are  discussed.  A  discussion  is  given  in  Annex  C  of  the  effect 
of  varying  some  of  the  assumptions  made  in  these  sections.  Results  are 
presented  of  the  ’gain"  obtained  when  different  options  are  used.  This 
‘gain*  is  defined  as  the  relative  improvement  in  P  over  its  value  in  the 
base  case  which  is  defined  in  the  next  paragraph.  This  base  case  value  is 
the  ratio  of  average  signal  power  to  average  noise  power  over  the  analysis 
penoa. 


The  base  case  arises  if  the  following  simple  (but  not  necessarily 
accurate)  model  of  reverberation  is  used.  Suppose  that 

E[  g(t.;\f.e*)  gfuJT.r.e-)  ]  -  (*\K2)  8<f-u)  5<r-p)  [4  4] 

so  that  all  the  reverberation  return  is  at  frequency  p,  reverberation  at 
different  phases  is  independent  and  the  reverberation  from  different 
scatter  ere  giving  contributions  at  the  same  phase  are  totally  correlated. 

In  this  case  for  the  unshaped  CW  pulse  (i.e.  W(t)-1 .  M(t)-O) 

E[  x(t)x(u)  ]  -  R  cos  2xp(t-u)  [4  5] 

and  the  scaling  factors  in  equation  [4.4]  are  chosen  so  that  the  average 
reverberation  power  is  R  .  If  no  shading  or  frequency  modulation  is  used 
(i.e  H(t)-l.N<t)-0), 

a2  -  R  ^  (T-v)  cos  2xfv  cos  2xpv  dv  [4  6] 

-B  f itfiit-ulT  +  ttfiltmlT  j  [4  7] 

4  l  X2  (f-u)2  X2  (f+u)2  J 

and 

Q1  -  ED  [4  8] 

Ignoring  terms  of  order  (2xfT)  \  which  is  assumed  small,  and  assuming 
that  the  reverberation  frequency  n  equals  the  analysis  frequency  f.  P  has 
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its  maximum  value  when  the  analysis  period  [O.T]  coincides  with  the 
entire  signal  period  [D,E]  so  that  D-0,  T-K  and  this  maximum  value  is 


which  is  the  average  signal  power  to  average  reverberation  power  ratio. 
This  is  the  base  case  with  which  the  other  options  are  compared  in  the 
remainder  of  this  report  in  order  to  assess  the  "gain*  achieved,  by  either 
doppler  processing  with  a  CW  pulse  or  using  an  FM  chirp  pulse  in  a 
reverberation  environment  with  given  characteristics.  The  average 
reverberation  power  R  is  proportional  to  the  pulse  length  K  (since  the 
width  of  the  scattering  region  is  cK/2,  where  c  is  the  speed  of  sound).  This 
must  be  taken  into  account  in  sections  4.1  or  4.2  or  Annex  C  when 
assessing  the  effect  of  a  changed  pulse  length. 

In  sections  4.1  and  4.2.  the  following  (more  realistic)  form  of  the 
covariance  of  g  is  used: 

E[  g(t.5’.r,8')  gfu.c’.r.e*)  ]-2R  L(f)  H f-n  «((;•-;•)  w-e*)  /  (*K) 

where  the  amplitudes  g  are  independent  for  different  phases,  centre 
frequencies  or  path  lengths  and  L(f)  gives  the  frequency  distribution  of 
the  reverberation  power  and  is  scaled  so  that  the  average  reverberation 
power  over  the  the  analysis  period  is  R.  Other  forms  will  be  considered  in 
Annex  C.  The  expected  value  of  x(t)x(u)  may  be  written  for  fcuit-K, 

E[  x(t)x(u)  ]-  (R/K)  M  L(f)  Win')  W(u-;'). 

cos(M(t-;,)-M(u-;')t-2itf  (t-u))  df  d£ 

For  u2Cu-K,  E[  x(t)x(u)  ]  is  given  by  an  equivalent  expression  with  u  and  t 
interchanged.  E[  x(t)x(u)  ]  is  zero  if  |t-u|>K.  In  the  cases  considered  in 
section  4.  it  is  shown  in  Annex  C  that 

E[  x(t)x(u)  ]  -  (R/K)  L(t-u)  W(t-u)  cos  2*n(t-u) 

where  the  functions  LO,  WO  and  the  quantity  u  are  defined  by  equations 
[C.24]  and  (C  .18}.  In  this  case. 
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a2  -  (R/K)  {  L(v)  W(v)  H(v)  cos  2*fv  cos  2*pv  dv  [4.1 3] 

-  (%K)  o  Mv)  W(v)  H(v)  cos  2jt(f-p)v  dv  [4.14] 

where  H()  is  defined  by  [C.40].  The  gain  G  may  therefore  be  written 
KQ,2 

G  -  T  [4.15] 

2  /  L(v)  W(v)  H(v)  cos  2*(f-p)v  dv 
o 

In  sections  4.1  and  4.2,  attention  will  be  confined  to  specific 
functions  M,  N,  L,  H  and  W.  The  latter  three  are  listed  in  tables  1  to  3.  Of 
the  functions  L  listed  in  table  1 ,  the  most  probable  to  occur  in  practice  are 
functions  1  or  4,  which  would  be  based  on  stationary  scatterers  or 
scatterers  whose  radial  components  of  velocity  relative  to  the  sonar  are 
normally  distributed.  For  each  of  the  functions  2  to  4,  the  parameter  a  is 
approximately  half  the  3dB  bandwidth  of  the  reverberation  power.  Each  of 
these  functions  approach  function  1  as  a-*0.  The  three  possible  functions 
W(t)  or  H(t)  are  denoted 

1)  NO  for  no  shading  or  shaping 

2)  EXP  for  exponential  (gaussian)  shading  or  shaping 

3)  COS  for  cos  squared  shading  or  shaping. 

Each  are  normalised  by  equations  [2.2]  or  [2.6]  and  each  are  expressed 
relative  to  a  time  origin  set  at  the  start  of  the  pulse  (W)  or  the  start  of 
the  analysis  period  (H).  The  functions  W  and  H  in  tables  2  and  3  are  valid 
for  CW  pulses  only,  where  M(t)-N(t)«0. 
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No 

1 

2 

3 

4 


Table  1 


L(f) 


The  Functions  L(f) 
L(v) 


S(Hi) 


1 


_ a 

Jt  ( (f-n)2  +  a2 ) 

&Z2^-Sin2gf-iiW0.72  ' 

jc  (f-(i)2 

Iln-2105fixp  (-( ff-u)2iin2]za2  )i 

wosa 


exp(-2nav) 

( 1  -  0.72nav )  if  |v|  <  (0.72:ta)'1 
0  if  | v|  >  (0.7271a)-1 

exp(  -Tc2  a2  v2  /(In  2)) 


No.  W(t) 

1  (NO)  1 


Table  2 

The  Functions  Wfti 

W(v) 


K-v 


2  (EXP)  Kexpt-rav/Kigiartnfi.v/ku 

n  ■  (erf(3))  3  erf(3) 

3  (COS)  a0  5ias2 izL(t/K-0,5) )  (K-v)  (2  +  cos  2itv/K)/3  + 

3  '  sin(27iv/K)/(27t) 
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Table  3 


No.  H(t) 

1  (NO)  1 


2  (EXP)  g0  5( 


Jt°- 25  (erf(3))0s 


3  (COS)  ao  sjfflS2in  tVT-0.5)  \ 
3o.s 


(T-v)  (2  +  cos  2jtv/T)/3  + 
sin(2jiv/T)/(2ji) 
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4.1  CW  Transmission 

If  the  transmitted  pulse  is  CW, 

M(t)  -0  [416] 

so  that 

W(v)  -  f  W(w-v)  W(w)  dw  [417] 

V 

t  K-v 

-  J  W,  (v-w)  W,  ( v+w)  dw  [4  18] 


where 

W^x)  -  W((K-x)/2)  [4.19] 

Assuming  that  N(t)-0,  which  gives  the  maximum  gain  against  noise,  a 
similar  expression  may  be  written  for  H(v)  with  H  replacing  W  and  T 
replacing  K  throughout.  The  expressions  for  H(v)  and  W(v)  from  tables  2 
and  3  may  then  be  used.  If  D-0  and  T-K, 

Q.  fTH(t)  W(t)  dt  [4  20] 

1  o 

and  the  equation  [4.15]  may  be  used  to  evaluate  the  gain.  This  gain  is 
presented  as  a  function  of  a  non-dimensional  parameter  co«27t(f-p)T  in 
figures  3  to  21  for  all  combinations  of  the  functions  L,  H  and  W  from 
tables  1  to  3.  Because  equation  [4.15]  is  symmetric  in  H  and  W  in  the  cases 
considered  here,  the  same  gain  is  obtained  if  the  shading  or  shaping 
functions  are  interchanged. 

The  following  features  may  be  noted  from  the  figures: 

(i)  In  figure  3,  which  presents  the  gain  in  the  case  that  all  the 
reverberation  is  confined  to  one  frequency  (function  1  from  table  1  for  L), 
there  is  some  gain  when  co-0  due  to  the  use  of  a  different  reverberation 
model  from  the  one  used  in  the  base  case.  The  difference  lies  in  the 
correlation  (or  otherwise)  of  contributions  to  the  reverberation  from 
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scatterers  with  different  path  lengths. 

(ii)  The  doppler  shift  of  the  signal  is  given  approximately  by 

Af  -2vF/c  [4.21] 

where  v  is  the  radial  component  of  the  relative  velocity  of  the  target  and 
sonar  and  c  is  the  speed  of  sound, 

to  -  (4jtFT/c)  v  [4.22] 

if  the  reverberation  is  centred  on  the  carrier  frequency  (ie  p=F)  and  the 
analysis  frequency  coincides  with  the  signal  (doppler  shifted)  frequency. 

Thus  a  and  hence  the  gain  (which  is  an  increasing  function  of  to  in  all 
figures)  will  be  larger  for  a  sonar  with  a  larger  carrier  frequency  F  if  the 
other  parameters  are  unchanged.  The  value  of  <o  will  also  be  larger  for  a 
longer  pulse,  but  the  effect  of  changing  T  on  the  gain  is  more  complicated 
for  two  reasons.  Firstly,  increasing  T  will  increase  the  average 
reverberation  power  R,  which  is  proportional  to  the  pulse  length  unless 
the  pulse  is  so  long  that  the  average  reverberation  power  from  equal  areas 
at  two  ranges  cT/2  apart  differs  significantly  from  one  another.  Thus  the 
detection  performance  will  only  improve  for  an  increased  pulse  length  if 
the  gain  improves  faster  than  linearly  with  to.  Secondly,  when  functions  2, 

3  or  4  are  used  for  L  (figures  4  to  21 ),  the  parameter  aT  is  proportional  to 
T,  so  that  increasing  T  will  move  from  one  line  to  another  giving  a 
possibly  reduced  gain. 

(iii)  Considerable  improvements  may  be  made  in  the  gain  for  a  fixed  value 
of  to  by  shaping  the  pulse  or  using  a  shading  function.  The  aim  of  shading  or 
shaping  is  to  reduce  the  discontinuities  in  the  reverberation  due  to  the 
sudden  turning  off  or  on  of  contributions  from  a  scatterer  when  t»£+K  or 

t»£  respectively  (shaping)  and  the  sudden  start  and  end  of  the  analysis 
period  (shading).  These  discontinuities  widen  the  reverberation  spectrum 
or  the  filter  spectrum  and  lead  to  a  reduced  gain.  The  function  1  (no 
shading  or  shaping)  has  large  discontinuities  at  the  ends  of  the  interval 
and  will  give  the  worst  gains.  The  function  2  (exp  shading  or  shaping)  has 
very  small  discontinuities  at  the  endpoints  and  gives  improved  gains.  The 
function  3  (cos  shading  or  shaping)  is  continuous  and  has  a  continuous 
derivative  at  the  endpoints.  The  discontinuity  in  the  second  derivative  at 
the  endpoints  yields  some  widening  of  the  spectra,  but  the  effect  is  much 
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smaller  than  in  the  other  cases.  Note  that  K  and  W  should  be  matched  to 
have  the  same  degree  of  continuity  as  there  is  little  to  be  gained  in 
improving  H  (say)  if  W  is  the  dominant  cause  of  the  decrease  in  gain.  The 
possibility  that  the  nature  of  the  propagating  medium  and  the  scatterers 
may  limit  the  reduction  in  the  discontinuities  and  therefore  the  gain 
whatever  the  functions  H  and  W  is  discussed  briefly  in  Annex  C. 

(iv)  The  effect  of  broadening  the  reverberation  spectrum  by  using 
functions  2,  3  or  4  from  table  1  (due  to  motion  of  the  scatterer  or 
uncompensated  motion  of  the  sonar)  depends  very  much  on  the  function 
selected.  The  reasons  for  this  are  discussed  in  section  5.  If  the  half 
bandwidth  a  is  increased,  the  gain  is  also  increased  when  co=0  due  to  the 
decrease  in  maximum  reverberation  power  brought  about  by  the  broadening. 
However  the  increase  in  gain  with  co  is  usually  slower  with  the  larger 
values  of  a,  due  to  the  increase  in  reverberation  power  in  the  main 
analysis  frequency  band.  This  slower  increase  is  more  marked  using 
function  2  or  3,  where  the  fall  of  power  with  frequency  away  from  the 
centre  frequency  is  slow,  than  using  function  4,  where  the  fall  is  much 
faster  beyond  the  3  dB  bandwidth.  In  fact  with  function  4  the  gain 
approaches  that  with  a=0  (function  1)  for  sufficiently  large  to,  indicating 
that  there  is  no  longer  any  appreciable  reverberation  power  in  the  main 
analysis  frequency  band. 
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Figure  17(11)  Gain  using 
Function  4  for  L 
Exp  Shading  and  Exp  Shaping 
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Figure  18(U)  Gain  using 
F unct ion  4  for  L 
Cos  Shading  and  Cos  Shaping 
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Figure  1S(U)  Gain  using 
Function  4  for  L 
Exp  Shading  and  No  Shaping 
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4.2  Chirp  FM  Transmission 

If  the  transmitted  pulse  is  a  chirp, 

M(t)  -  2*p(t-K/2)2  OStSK  [4.23] 

where  2{5K  is  the  FM  bandwidth,  so  that 
.k 

W(v)  -  S  W(w-v)  W(w)  cos  (  2rc|}(  2vw  -  v(v+K) ) )  dw  [4.24] 


-  J  W,(v-w)  W,(v+w)  cos  2jt|ivw  dw 


[4.25] 


where 

W,(x)  -W((K+x)/2)  -  W((K-x)/2)  [4.26] 

Assuming  that  N(t)-cos  2jt(l(t-T/2)2,  which  gives  the  maximum 
contribution  from  the  target  echo,  a  similar  expression  may  be  written  for 
H(v)  with  H  replacing  W  and  T  replacing  K  throughout.  However,  the 
expressions  for  H(v)  and  W(v)  from  tables  2  and  3  are  no  longer  valid.  If 
D-0  and  T-K, 

Q,  JoTH(t)W(t)dt  [4.27] 

and  the  equation  [4.15]  may  be  used  to  evaluate  the  gain.  This  gain  is 
presented  as  a  function  of  y  -  2n$Tz  (which  is  n  times  the 
bandwidth-pulse  length  product)  in  figure  22  for  all  combinations  of  the 
functions  H  and  W  from  tables  2  and  3  and  function  1  for  L  from  table  1 . 

There  is  very  little  difference  in  the  gain  for  the  other  functions  L  except 
for  small  values  of  y.  Again,  interchanging  the  shading  and  shaping 
functions  has  no  effect  on  the  gain. 

From  figure  22,  the  gain  is  an  increasing  function  of  y.  The  gain  is 
proportional  to  y  at  all  but  very  small  values  of  y  as  evidenced  by  the  3  dB 
increase  in  gain  (in  decibels)  for  each  doubling  of  y.  Thus  the  gain  is 
proportional  to  the  FM  bandwidth  and  to  the  pulse  length.  However,  as 
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noted  in  section  4.1 ,  the  reverberation  power  will  normally  be  proportional 
to  the  pulse  length  and  the  detection  performance  will  therefore  not  be 
improved  by  increasing  the  pulse  length. 
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5.  Approximate  Expressions  for  the  Gain 
5.1  CW  Pulses  with  large  to 

In  the  cases  considered  in  section  4.1 ,  the  gain  G  may  be  written 

G  -  <3, 2  /  2  ^  S(v)  cos  aw  dv  [5.1] 

where 


S(v)  -  L(vT)  W<vT)  H(vT). 


[5.2] 


For  a>  sufficiently  large  and  assuming  that  S  is  infinitely  differentiable, 
the  integral  in  the  denominator  in  equation  [5.1]  (which  will  be  denoted  J) 
may  be  expressed,  after  repeated  integration  by  parts, 


J 


f  S(2n)(1)  sin  <o  s<2n+1>(1)  cos  a)  -  S<2n+1)(0) 
X(-l)n  -  +  - 

n-C  ai2"4'1  C02n+2 


j  [5.3] 


where  S(t<)  denotes  the  kth  derivative  of  S.  J  will  be  dominated  by  the  first 
term  in  the  series  for  which  the  numerator  is  non-zero.  From  equations 
[C.18]  and  [C.34],  W(T)  -  H(T)  -  0,  so  that  S(1)  -  S’(1)  -  0.  If  S'(0)  *  0  (it 
must  be  negative  since  S(x)  is  decreasing  when  x  >  0), 

J  -  -S'fOJ/oj2  -  -T2  (  T  L'(0)  -t-  W(0)  +  H'(0) )/  co2  [5.4] 

for  large  <u.  After  differentiating  equation  [C.18]  and  noting  that  W  and  H 
are  symmetric  about  t  -  T/2, 

W(0)  -  -W2(0).  [5.5] 

Similarly, 

H*(0)  -  -H2(0).  [5.6] 

Consider  now  each  of  the  functions  L  listed  in  Table  -1 .  If  function  1  is 
used,  all  the  derivatives  of  L  are  zero  and 


42 


UNCLASSIFIED 


UNCLASSIFIED 


J  -  T2  (W2(0)+H2(0))/oj2  [5.7] 

provided  that  the  term  on  the  right  hand  side  is  non-zero.  This  will  be  true 
if  there  is  either  no  or  exp  shading  or  shaping  but  not  if  there  is  both  cos 
shaping  and  shading.  Notice  that  in  Figure  3,  the  gain  increases  by  6  dB  per 
doubling  of  co  indicating  that  G  is  proportional  to  co2  for  large  co  except  in 
the  cos  shading  and  shaping  case.  If  there  is  cos  shading  and  shaping, 

J  -  T6  (W"2(0)+H"2(0))/co8  [5.8] 

and  the  gain  increases  by  18  dB  per  doubling  of  co. 

If  function  4  is  used,  all  odd  derivatives  of  L(u)  are  zero  when  u-0  and 
all  derivatives  of  L(u)  are  small  when  u-T.  For  sufficiently  large  co,  the 
gain  will  approach  the  value  obtained  above  in  the  function  1  case; 
however  co  will  now  have  to  be  larger  before  this  limit  is  achieved  because 
the  even  derivatives  of  L(u)  at  u-0  can  be  quite  large,  especially  for  large 
values  of  a.  These  large  values  increase  the  odd  derivatives  S(2n+1)(0)  for 
nz2  so  that  a  larger  value  of  co  is  required  before  the  term  S(3)  (0)/co6  (say) 
may  be  neglected.  This  delay  in  reaching  the  limit  may  be  seen  in  figures 
1 6  to  21 .  For  large  a  and  moderate  co,  with  no  shading  or  shaping, 


J  -  T2  exp  (-pv2)  cos  cov  dv 
-  T2/2  V(7t/p)  exp  (-o^Mp) 

where  p  -  it2  a2  T 2  /  (In  2),  which  is  quite  large.  Hence  the  gain  is 
proportional  to  exp  (co 2/4p),  which  gives  the  parabolic  shape  of  the  gain  in 
decibels  for  the  smaller  values  of  co  in  figures  16  to  21 .  Note  that  when  co 
is  not  too  large,  the  gain  is  reasonably  independent  of  the  shading  or 
shaping  functions  used. 

If  function  2  is  used, 

L'(0)  -  -2rta  [5.9] 

Hence  S'(0)*0  whatever  shading  or  shaping  functions  are  used  and  |S'(0)| 
will  be  larger  and  hence  the  gain  less  as  a  increases.  The  gain  will  still 
increase  at  the  rate  of  6  dB  per  doubling  of  co.  This  may  be  seen  in  figures 
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4  to  9. 

If  function  3  is  used,  equation  [5.3]  may  not  be  used  as  S  is  not 
differentiable  when  1.18  avT  -  1,  which  occurs  within  the  integration 
interval  in  the  cases  considered  in  section  4.1 .  The  upper  limit  for  the 
integral  in  equation  [5.1],  1,  must  be  replaced  by  (0.72  naT)'1 ,  which  will 
be  denoted  x-  Thus, 

J  -  (S’(x)  cos  (cox)  ■  S'{0))  /co2 

-  (L’(x)  H(x)  W(x)  cos  (ox)  -  L*(0)  -  W(0)  -  H'(0) )  /a>2 

This  has  an  oscillating  term  (with  decreasing  frequency  but  increasing 
amplitude  for  increasing  a)  which  is  non-zero  for  any  of  the  shading  or 
shaping  functions  as  may  be  seen  in  figures  10  to  15.  The  upper  and  lower 
limit  of  the  oscillation  grow  with  a>2. 


44 


[5.10] 

[5.11] 


UNCLASSIFIED 


UNCLASSIFIED 


5.2  FM  Pulses  with  large  v 


If  no  shading  or  shaping  is  used  and  the  reverberation  is  of  the  form 
given  by  function  1  in  table  1  with  f  -  p  then 


L(vT)  -  1 

[5.12] 

W(vT)  -  T  sin  (yv(l-v))  /  (yv) 

[5.13] 

H(vT)  «  T  sin  (yv(l-v))  /  (yv) 

[5.14] 

and 

J  -  T2^0sin2(Yv(1-v))/(yv)2dv 

[5.15] 

»  T2  y’1  0  sin2(u(1  -u/y))  /  u2  du 

[5.16] 

on  letting  u  «  yv_  Now  the  integrand  in  equation  [5.16]  is  bounded  in 
magnitude  by  an  integrable  function  (1  for  0£u£l  and  1/u2  for  u>1 )  and 
approaches  sin2  u  /  u2  as  y-*».  By  the  Lesbegue  dominated  convergence 
theorem  (see  reference  [2],  page  53) 

lim  rYsiniu(1-u/y)  du 

J0  u2 

lim  sinLuCl-u/y)  du 
o  y-*°  u2 

[5.17] 

'j 

r°°  sinludu 
lo  u2 

[5.18] 

-  n/2. 

[5.19] 

Thus 

J  *  3tT2/2y+0(y15) 

[5.20] 

The  next  term  in  the  asymptotic  series  is  Ofy'1,5)  since  on  letting  u  -  yV 
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lim  y< 

y-*oo 


/. 


YsiQiug^i^siniudu 


[5.21] 


_  lim 

y— k» 


A  1lnyiUj 

Jo  u2 


This  limiting  behaviour 
If  other  functions 

since 


»  0  if  k  <  0.5. 
may  be  seen  in  figure  22. 

L  are  used,  the  limiting  behaviour  is  the  same, 


[5.22] 


[5.23] 


r' 

J  "  T  o  L(u/y)  sin2(u(1-u/y))  /  u2  du 


25?'  ~  * 


integrand  as  y->~  /s  still  rc/2.  Equation  [5.20]  still  applies. 

If  more  general  shading  or  shaping  functions  are  used, 


and  the  limit  of  the 


[5.24] 


L(vT/y)  W(vT/y)  H(vT/y)  dv 


so  that 

J  ~  x  T2/y 


where 


[5.25] 


[5.261 


X 


'  T'2  JQ  L(vT/y)  W(vT/y)  H(vT/y)  dv. 

r-*~ 


[5-27] 


Now 


46 

UNCLASSIFIED 

f 


\ 


UNCLASSIFIED 


T-v/y 


W(vT/y)  =  T W^Tv/y-Tw)  W^Tv/y+Tw)  COS  vw  dw 
r1 

T  0  W12(Tw)  cos  vwdw  =  h,(v)  (say) 


as  y-*».  Similarly, 


[5.28] 

[5.29] 


H(vT/y)  -  T  | 


i-v.y 


->  T 


r 


H^Tv/y-Tw)  H,(Tv/y+Tw)  cos  vw  dw 
H12(Tw)  cos  vw  dw  =  h2(v)  (say) 


[5.30] 

[5.31] 


as  y->~.  Since  |W(vT/y)|  and  |H(vT/y)|  are  each  bounded  by  a  constant  times 
v  for  large  v,  the  product  |W(vT/y)  H(vT/y)|  is  bounded  by  an  integrable 
function  and  the  Lesbegue  dominated  convergence  theorem  implies  that 


X  =  S  h^v)  h2(v)  dv. 
r1 

=  (Tt/2)  J0  H^Tw)  W12(Tw)  dw 

The  equality  of  the  integrals  in  equations  [5.32]  and  [5.33]  is  demonstrated 
in  annex  C.3,  where  it  is  also  shown  that  the  gain  has  its  maximum  value 
for  a  large  fixed  y  when  no  shading  or  shaping  are  used.  The  gain  is 
proportional  to  y  and  therefore  increases  linearly  with  bandwidth. 
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6.  An  Example 

To  illustrate  the  use  of  the  results  presented  in  this  paper,  a  simple 
example  of  a  sonar  attempting  to  detect  a  target  is  considered.  In  this 
section,  the  average  signal  power  P,  the  background  noise  in  a  1Hz  band,  N0 
and  the  average  reverberation  power  R  will  be  expressed  in  decibels  (re 
Ippascal)  and  may  be  calculated  using  the  usual  sonar  equations  as 

P  -  SL-2PL  +  TS 

N0  -  NL  -  Dl 

R  -  SL  -  2  PL  +  S  +  1 0  log  A 

where  SL  is  the  source  level,  PL  is  the  one  way  propagation  loss  to  the 
target,  TS  is  the  target  strength,  NL  the  ambient  noise  level,  Dl  the 
directivity  index,  S  the  scattering  strength  and  A  the  reverberating  area. 

In  this  example,  the  sonar  transmits  2  second  pulses  at  4  kHz  and  the 
threshold  is  set  to  give  a  false  alarm  probability  of  1 0'5.  Suppose  that  for 
the  analysis  period  during  which  the  target  signal  is  maximum, 

P  *  60  dB 

N0  -  48  dB 

R  -  85  dB 

Consider  first  the  noise  limited  performance  of  the  sonar  ie  the 
performance  of  the  sonar  in  the  presence  of  noise  only.  If  the  conditions  of 
section  3  to  give  a  maximum  value  of  P  are  met, 

P  -  PT/2N0  -  12  dB.  [6.7] 

Hence,  using  equation  [2.8]  or  figure  1 ,  the  probability  of  detection  is 
about  0.8.  The  signal  excess,  which  is  defined  as  the  ratio  of  P  to  the  value 
of  P  required  for  a  50%  probability  of  detection  (10  dB  here)  is  about  2  dB. 

Next  consider  the  reverberation  limited  performance  of  the  sonar  ie 
the  performance  of  the  sonar  in  the  presence  of  reverberation  alone.  If  the 
conditions  of  the  base  case  of  section  4  apply, 

P  -  P/R  -  -25  dB  [6.8] 


[6.4] 

[6.5] 

[6.6] 


[6.1] 

[6.2] 

[6.3] 


48 


UNCLASSIFIED 


UNCLASSIFIED 


and  the  probability  of  detection  is  1 .04x1  O'5,  only  marginally  above  the 
false  alarm  probability.  The  signal  excess  is  about  -35  dB.  If  doppler 
processing  ,  as  described  in  section  4,  is  used,  table  4  gives  the  signal 
excess  and  probability  of  detection  for  three  combinations  of  shading  and 
shaping  functions  and  several  values  of  the  radial  component  of  the 
velocity  of  the  target  relative  to  the  sonar  on  the  assumption  of  no 
reverberation  frequency  spreading.  These  are  derived  by  first  calculating 
the  doppler  shift  Af  induced  in  the  signal  by  the  radial  speed  ( -  2.8  x 
radial  speed  (in  knots)  for  4  kHz  sonar)  and  using  this  to  calculate  co  = 

2rcAfT  and  hence,  using  figure  3,  the  gain  in  P  over  its  base  case  value.  The 
signal  excess  is  this  improved  value  (in  decibels)  less  10  dB  as  required 
for  the  50%  probability  of  detection.  The  probability  of  detection  is  found 
from  figure  1 .  Note  that  even  with  zero  dopploer,  there  is  a  slight 
improvement  in  signal  excess  over  the  base  case  value  due  to  the  different 
assumptions  made  about  the  nature  of  reverberation. 

Against  a  combined  noise  and  reverberation  background,  the  high 
values  of  signal  excess  and  detection  probabilities  above  0.8  would  not  be 
achievable.  Instead  a  value  of  P  given  in  terms  of  the  value  for  noise 
limited  performance  Pn  and  the  value  for  reverberation  limited 
performance  Pr  by 

P-1  -  Pn  i  +  Pr-i  [6.9] 

should  be  used  to  calculate  the  signal  excess  and  probability  of  detection. 

Table  5  gives  the  signal  excesses  and  the  probabilities  of  detection 
against  the  combined  background.  Note  that  the  signal  excess  can  never  be 
better  than  the  value  calculated  against  noise  alone  or  the  value 
calculated  against  reverberation  alone. 

If  there  is  a  spread  of  reverberation  frequencies,  then  figure  3  is  not 
the  appropriate  one  to  calculate  the  gain.  If,  for  example,  the  spread  in  the 
reverberation  frequencies  had  a  bandwidth  of  4  Hz  and  a  gaussian  shape, 
then  function  4  from  table  1  with  aT-4  (a  is  the  half-bandwidth)  should  be 
used.  The  signal  excesses  and  the  probabilities  of  detection  in  this  case, 
using  the  aT-4  lines  in  figures  16  to  18  to  calculate  the  gain,  are  given  in 
table  6.  Note  that  there  is  no  longer  any  performance  with  only  one  or  two 
knots  of  relative  radial  speed  with  exp  or  cos  shaping  and  shading. 

If  an  FM  chirp  pulse  with  bandwidths  200,400  or  800  Hz  is  used, 
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instead  of  the  CW  pulse  with  doppler  processing,  the  gain  may  be 
calculated  from  figure  22  with  7  -nxbandwidthxpulse  length.  The  signal 
excesses  and  probabilities  of  detection  are  given  in  table  7  and  these  are 
not  sensitive  to  frequency  spreading  of  the  reverberation.  Note  that  this 
processing  option  does  not  give  any  useful  performance  in  this  case. 

Table  4 

Signal  Excesses  (dB)  and  Probabilities  of  Detection  with  Reverberation 
Limitations  Only  (CW  Pulse) 


Relative  Radial 
Speed  (knots) 

No  Shading  or 
Shaping _ _ 

Exp  Shading  and 
Shaping _ 

Cos  Shading  and 
Shaping 

SE 

- - - 

Dd 

SE 

pd 

SE_  _ 

pd 

0.0 

-33 

0.00001 

-31 

0.00001 

-32 

0.00001 

0.5 

-16 

0.00006 

-13 

0.0002 

-3 

0.07 

1.0 

-10 

0.0007 

22 

1. 

19 

1. 

2.0 

-4 

0.03 

30 

1. 

39 

1. 

4.0 

2 

o.a 

36 

1. 

57 

1. 

Table  5 


Signal  Excesses  (dB)  and  Probabilities  of  Detection  with  Reverberation 
and  Noise  with  no  Frequency  Spreading  (CW  Pulse) 


Relative  Radial 
Speed  (knots) 

No  Shading  or 
Shaping _ 

Exp  Shading  and 
Shaping 

Ccs  Shading  and 
Shaping 

SE 

Dd 

.  SE 

od 

SE 

Dd 

0.0 

-33 

0.00001 

-31 

0.00001 

-32 

0.00001 

0.5 

-16 

0.00006 

-13 

0.0002 

-4 

0.07 

1.0 

-10 

0.0007 

2 

0.8 

2 

0.8 

2.0 

-5 

0.02 

2 

0.8 

2 

0.8 

4.0 

-1 

0.2 

2 

0.8 

2 

0.8 
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Table  6 

Signal  Excesses  (dB)  and  Probabilities  of  Detection  with  Reverberation 
and  Noise  and  Frequency  Spreading  (CW  Pulse) 


Relative  Radial 

No  Shading  or 

Exp  Shading  and 

Cos  Shading  and 

Sneed  fknotsi 

Shaoino 

Shaoino _ 

Shaoino _ 

SE 

od 

SE 

Dd 

SE 

od 

0.0 

-26 

0.00001 

-26 

0.00001 

-26 

0.00001 

0.5 

-24 

0.00002 

-24 

0.00002 

-24 

0.00002 

1.0 

-20 

0.00002 

-20 

0.00002 

-20 

0.00002 

2.0 

-8 

0.002 

-5 

0.02 

-5 

0.02 

4.0 

-1 

0.2 

2 

0.8 

2 

0.8 

Table  7 


Signal  Excesses  (dB)  and  Probabilities  of  Detection  with  Reverberation 
and  Noise  (FM  pulse) 


Bandwidth  (Hz) 

No  Shading  or 
Shaoino 

Exp  Shading  and 
ShaDina 

Cos  Shading  and  1 
ShaDina 

SE 

od 

SE 

Dd 

SE 

od 

200 

-9 

0.001 

-13 

0.0002 

-12 

0.0003 

400 

-7 

0.004 

-10 

0.0007 

-9 

0.001 

800 

-4 

0.03 

-8 

0.002 

-7 

0.004 
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LSummaa 

The  probability  of  detection  of  a  reflected  signal  from  a  target  in  the 
presence  of  interference  depends  only  on  the  quantity  P,  defined  by 
equation  [2.9],  for  a  given  probability  of  false  alarm  provided  certain 
assumptions  listed  in  annex  B  regarding  the  nature  of  the  interference  are 
met  This  relationship  is  illustrated  in  figure  1 .  Most  of  the  variations  in 
the  assumptions  that  are  considered  in  annex  B  and  illustrated  in  figures 
23  to  28  do  not  change  the  value  of  P  required  to  achieve  a  particular 
probabilty  of  detection  by  more  than  a  few  decibels.  The  variation  with  the 
greatest  effect  is  the  rejection  of  the  gaussian  assumption. 

The  quantity  P  may  be  estimated  if  the  statistics  and  other 
properties  of  the  interference  are  known.  However,  P  is  highly  dependent 
on  the  following  factors 

i)  .  the  pulse  shaping  and  analysis  filter  shading  functions  used 

ii)  the  statistics  of  the  reverberation 

iii)  the  frequency  spread  of  the  reverberation  and  the  difference  between 
the  main  reverberation  and  the  signal  centre  frequencies  if  doppler 
processing  is  used 

iv)  the  bandwidth  if  FM  pulses  are  used 

v)  the  alignment  of  the  pulse  and  analysis  periods  and  the  alignment  of 
the  pulse  and  analysis  frequencies 

vi)  the  relative  values  of  the  pulse  length  and  the  length  of  the  analysis 
period. 

In  the  presence  of  white  gaussian  noise,  the  maximum  value  of  P  is 
given  by 

P  »  PT/2Nq 

where  P  is  the  average  signal  power,  T  is  the  pulse  length  and  N0  is  the 

noise  power  spectrum  level.  This  leads  to  the  familiar  expressions  for  the 
detection  threshold  (eg  reference  [1]  Chapter  12) 

Using  specific  options  regarding  the  factors  listed  above,  a  base  case 
for  detection  in  the  presence  of  reverberation  may  be  defined  where  P  is 
simply  the  ratio  of  the  average  signal  power  to  the  average  reverberation 
power.  This  base  case  leads  to  very  poor  performance  in  most  cases,  but 
this  performance  may  be  improved  by  using  other  options.  The  resulting 
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proportional  increase  in  P  is  denoted  the  gain. 

Assuming  a  particular  model  of  the  statistics  of  the  reverberation,  an 
expression  for  the  gain  is  derived  (equation  [4.1 5])  as  a  function  of  the 
other  factors  listed  above  as  i)  and  iii)  to  vi). 

This  expression  is  used  to  calculate  the  gain,  assuming  no 
misalignment  (see  v)  above)  and  equal  pulse  and  analysis  periods  using 
three  options  each  of  shading  and  shaping  functions.  For  doppler  processing 
in  the  case  of  no  reverberation  frequency  spreading,  the  gain  is  presented 
in  figure  3  as  a  function  of  co  «  2nAfT,  where  Af  is  the  difference  between 
the  signal  frequency  and  the  reverberation  frequency  and  T  is  the  pulse 
length.  Similar  results  are  presented  in  figures  4  to  21  for  various  options 
of  reverberation  frequency  spreading  and  illustrate  the  importance  of  this 
spreading  and  the  large  differences  in  the  results  which  may  be  obtained  if 
different  assumptions  are  made  about  the  shape  of  the  spreading  function. 

The  spreading  could  be  caused  by  motion  of  the  scatterer  relative  to  the 
sonar  (eg  sea  surface  reverberation)  or  if  the  sonar  is  moving  by  the  fact 
that  all  reverberation  entering  a  beam  will  not  have  the  same  doppler 
shift.  The  same  effect  could  be  caused  if  the  scatterers  change  their 
properties  over  time  periods  less  than  a  pulse  length.  Figure  22  shows  the 
gain  using  FM  chirp  pulses  as  a  function  of  y-  JtWT  where  W  is  the 
bandwidth.  This  gain  is  not  sensitive  to  the  reverberation  frequency 
spreading  considered  in  the  doppler  processing  case.  The  gain  is  greatest 
when  no  shading  or  shaping  is  used. 

Asymptotic  expressions  for  the  gain  for  large  values  of  co  using 
doppler  processing  or  large  values  of  y  using  FM  pulses  are  derived  in 
section  5.  It  is  shown  that  under  most  circumstances  the  gain  is 
proportional  to  co2  for  large  co  using  doppler  processing  (ie  6  dB 
improvement  in  gain  for  each  doubling  of  co)  and  the  gain  is  proportional  to 
yfor  large  y  using  FM  pulses.  The  proportionality  constant  depends  on  the 
shading  and  shaping  functions  used  and  for  doppler  processing  on  the 
function  assumed  for  the  frequency  spreading  of  the  reverberation. 

The  combined  effect  of  noise  and  reverberation  is  to  give  a  value  of  P 
calculated  as 

P'1  .  Pn-i  +  P-'  [7.2] 

where  Pn  is  the  value  of  P  calculated  with  interference  due  to  noise  only 


53 


UNCLASSIFIED 


UNCLASSIFIED 


and  Pris  the  value  of  P  calculated  with  interference  due  to  reverberation 
only. 


The  mis-alignment  of  the  pulse  and  analysis  periods  for  doppler 
processing  can  give  quite  large  degradations  (up  to  20  dB  using  exp  shading 
and  shaping  as  defined  in  tables  2  and  3)  in  the  gain.  This  degradation  may 
be  reduced  by  overlapping  the  analysis  periods,  which  necessitates  more 
processing  and  may  increase  the  false  alarm  rate.  Misalignment  of  the 
pulse  and  analysis  frequencies  using  doppler  processing  produces  a 
smaller  degradation  (at  worst  4  dB  using  no  shading  or  shaping)  assuming 
that  the  analysis  frequencies  are  separated  by  T1. 

The  misalignment  of  the  pulse  and  analysis  periods  for  FM  pulses  is 
analogous  to  misalignment  of  the  pulse  and  analysis  frequencies  for 
doppler  processing.  Providing  the  analysis  periods  are  overlapped  so  that 
their  starting  times  are  separated  by  W'1,  the  maximum  degradation  is  4 
dB  when  no  shading  or  shaping  is  used.  The  misalignment  of  the  pulse  and 
analysis  frequencies  for  FM  pulses  is  analogous  to  misalignment  of  the 
pulse  and  analysis  periods  for  doppler  processing.  Although  the  selection 
of  a  different  analysis  period  may  reduce  the  frequency  misalignment,  the 
degradation  may  be  quite  large  (up  to  20  dB  for  exp  shading  and  shaping  as 
defined  in  tables  2  and  3). 
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ANNEX  A 

PROBABILITY  OF  FALSE  ALARM  AND  PROBABILITY  OF  DETECTION 

When  a  signal  is  not  present,  the  output  y(t)  from  the  beam-former  is 
made  up  entirely  of  interference  x(t).  Assume  that,  in  this  case,  B  and  C, 
as  defined  in  equations  [2.4]  and  [2.5],  have  a  Gaussian  distribution  with 
zero  mean.  The  variance  of  B  is  then 

var(B)  mf  J  S(t,u)  cos  (N(t)+2rcft)  cos  (N(u)+2itfu)  dt  du  [A.1] 

o  o 

where 

S(t,u)  -E[x(t)x(u)]H(t)H(u).  [A.2] 

On  letting  v-t-u  in  equation  [A.1], 

var(B)  -  /  \  S(t,t-v)  { cos  (N(t)-N(t-v)+27ifv) 

+  cos  (N(t)+N(t-v)+2nf(2t-v)) }  dt  dv.  [A.3] 


Similarly, 

T  T 

var(C)  -  i  J  S(t,t-v)  { cos  (N(t)-N(t-v)+2rtfv) 

-  COS  (N(t)+N(t-v)+27tf(2t-v)) }  dt  dv  [A.4] 

.T  T 

cov(B.C)  -  Jq  S(t.t-v)  sin  (N(t)+N(t-v)+2jrf(2t-v))  dt  dv.  [A.5] 

On  making  the  substitution  w  -  T  -  2t  +  v  in  the  second  of  the  inner 
integrals  in  equation  [A.3]  or  [A.4], 

.t 

I  S(t,t-v)  cos  (N(t)+N(t-v)+27tf(2t-v))  dt  - 

v 

.  T-v 

0.5  J  S,(v-w,v+w)  cos  (N1(w-v)+N1(w+v)+27tf(T-w))  dw  [A.6] 

where 
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S^x.y)  -  S(  (T+x)/2  ,  (T-y)/2 ) 

[A.7] 

N^x)  -  N((T-x)/2). 

fA.8} 

Thus, 

1 

,T  T 

J  /  S(t,t-v)  cos  (N(t)+N(t-v)+  2jrf(2t-v))  dt  dv 

0  v 

■  a,  cos  2jcfT  +  a 2  sin  2rcfT 
»  a  cos  (2jcfT  -  £j) 

[A.9] 

where 

T  T-v 

ai  *  0.5  Jq  \  S,  (v-w,v+w)  cos  (27cfw-N1  (w-vJ-N^w+v))  dw  dv 
-  a  cos  Cj 

[A.  10] 

T  T-v 

“0-5  M^S^v-w.v+w)  sin  (2jtfw-N1(w-v)-N1(w+v))  dwdv 
«  a  sin  Ej 

[A.  11] 

Similarly, 

fT,  T 

i  /  S(t,t-v)  sin  (N(t)+N(t-v}+  2icf(2t-v))  dt  dv 

0  v 

-  a,  sin  2«fT  -  ag  cos  27ifT 

-  a  sin  (27tfT  -  Ej) 

[A- 12] 

Thus, 

var  (B)  -  o2  (1  +  A  cos  e) 

[A.13] 

var  (C)  -  o2  (1  -  A  cos  e) 

cov(B.C)  »  o2  A  sin  e, 

where 

[A.  14] 

[A.  15] 

r 

tV 

o2  -  J0  Jv  S(t,t-V)  cos  (N(t)-N(t-v)+2nfv)  dt  dv 

{A.  16] 
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J  T-v 


■  0.5  J  /  S,(v-w,v+w)  cos  (2jrfv+N1(w-v)-N1(w+v))  dw  dv 

0  -T+v 


V 


e  «  2jtfT  -  e 


The  correlation  coefficient  between  B  and  C  is 


A  sine 


(1  -  A2  cos2  e) 0  5 

and  the  joint  probability  density  function  for  B  and  C  is 


[A.  17] 
[A- 18] 
[A- 19] 

[A.20] 


f(B,C)  « 


_L 


_  oxp  (  B2  +  C2  -  Af  (B2-C21  cos  e  ±  2BC  sin  el } 

2ji(1-A2)°'5o2  2a2  (1-A2) 


On  transforming  to  polar  co-ordinates 
r  «(B2  +  C2)05 
0  -  arctan(  C/B) 

the  joint  probability  density  function  of  r  and  0  becomes 

f(r,e) - 1 - exp  ( -  dll  -Acoa-(2fcel } . 

2ji(1-A2)0  5o2  2  o2  (1-A2) 

On  integration  with  respect  to  0,  the  probability  density  function  for  r 
becomes 


[A.21] 

[A.22] 

[A.23] 

[A.  24] 


f(r)  -■ 


_L 


o2  (1-A2)0  5 


exp(  — 


■)  I  (• 


JuL 


2  c2  (1-A2)  0  2  o2  (1-A2) 


)• 


[A.25] 


where  l0( )  is  the  modified  Bessel  function  of  order  zero,  defined  for 
example  in  reference  [3]. 
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If  |A|  ■  1 ,  the  above  analysis  breaks  down.  In  this  case,  the  correlation 
coefficient  between  B  and  C  is  1  or  -1  depending  on  the  sign  of  sin  e  and 

var(B)  -  2  o2  cos2  (e/2)  [A.26] 

var(C)  -  2  o2  sin2  (e/2).  [A.27] 

Thus 

C  =  B  tan  (e/2)  [A. 28] 

r  -  (  B2  +  C2  )0-5  -  - LEI -  .  - LEI -  [a.29] 

I  cos  (e/2)|  |sin  (e/2)| 

and  the  probability  density  function  for  r  is 

f(r)  -  — 1 —  exp  ( - ) .  [A. 30] 

V(2jc)  a  4  o2 

In  any  case,  if  R  is  the  threshold  described  at  the  end  of  the  first 
paragraph  in  section  2,  the  probability  of  false  alarm, 

pta  -f  f(r)dr.  [A.31] 

If  A  «  0  (which  is  approximately  true  in  most  cases  considered  in 
subsequent  sections),  the  probability  density  function  is 

f(r)  -  exp  ( — tL)  [A.32] 

o2  2c2 

and  the  probability  of  false  alarm  is 
2 

Pfa  "  exp  (--ft—).  [A.33] 

202 
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Now  suppose  that  a  target 
present.  The  variances  of  B  and 
same  as  in  the  interference  only 
The  expected  values  of  B  and  C 


echo  of  the  form  given  in  equation  [2.7]  is 
C  and  the  covariance  of  B  and  C  will  be  the 
case  as  given  by  equations  [A. 4]  to  [A. 6], 
will  no  longer  be  zero,  but 


fE 

E[B]  -  (»)0.s  lo  W(t.D)  cos  (M(t.D)  +  2nfif  +Qi)  H({)  cqs  (N(t)+2jtft)  dt 

*  (Vs H(t)  { cos  (M(t-D)-N(t)+2jr(f1-f)t+ei)  + 

[A. 34] 

cos  ( M(t-D)+N(t)+2js(f,+f)t+01) }  dt 

"  (P/2)°'s  ( Q,cos  (0,  +  a)  +  Q2  cos(  0,  +  3 ) } 

[A.35] 

[A.36] 

where 

Ql  C0S  a  ’  fD  W(t-D)  H(t)  cos  (  M(t-D)  -  N(t)  +  2n(frf)t )  dt 

[A.37] 

°1  sin  a  mID  W(t'D)  H(t)  sin  ( M(t-D)  -  N(t)  +  2n(f,-f)t )  dt 

[A.38] 

Q2  cos  a  =.  io  W(t-D)  H(t)  cos  (  M(t-D)  +  N(t)  +  2ji(f1+f)t )  dt 

[A. 39] 

Q2  sin  a  -  fo  W(t-D)  H(t)  sin  ( M(t-D)  +  N(t)  +  2n(f, +f)t )  dt. 

[A.40] 

Similarly, 

E[C]  "  (V5  { -Qis'n  (e,  +  a)  +  Q2  sin  ( 8,  +  (3 ) } 


[A.41J 


The  joint  probability  density  function  of  B  and  C  with  sional  Dresent 


exp  { -  &Af((B-P  ]2-fC-e  „  +2lB-pr)ir.-n)^m  f| 

- - - 2^(1-a^  _ _ 

2k  a2  (1  -  A2}0  5 


[A.  42] 
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On  transforming  to  polar  co-ordinates  using  equations  [A. 22]  and  [A.23], 
the  various  terms  in  the  above  equation  may  be  written 

(B-B  )2  +  (C -Cf  *  r2  -  (2P)as  r  (Q,  cos  (0+6,  +a)  +  Q2  cos  (6-0,  -p)) 

+(P/2)  (Q,2  +  2  Q,Q2  cos  (20,+a+P)  +  Q22)  [A.43] 

(B-B)2-  (C-C)2  =  r2  cos  20  -  (2P)as  r  (  Q,  cos  (0-6, -a)  +  Q2cos  (0+0,  +p) ) 

+(P/2)  (Q,2cos  (20,  +2a)  +  2  Q,Q2  cos  (a-p)  +  Q22  cos  (20,+2p)) 

[A.  44] 

2  (B-B )  (C-C )  -  r2  sin  20  -  (2P)as  r  (  Q,  sin  (0-0, -a)  +  Q2  sin  (0+0,  +p) ) 

+(P/2)  (-Q,2sin  (20,  +2a)  -  2  Q,Q2  sin  (a-p)  +  Q22  sin  (20,+2p)) 

[A.45] 

arid  the  joint  probability  density  function  for  r  and  0  is 
r  exp  ( - ^ - ) 

f(r,0)  =  2  q2  (1  -  A2 ) _  [A.46] 

2rr  o2  (1  -  A2)0  5 


where 

XX  -  1^(1 -A  cos  (20-e))  -  (2P)0-5  r  (Q,  cos  (6+0, +a)  +  Q2cos  (6-6,-P)) 

+A  (2P)as  r  (Q,  cos  (0-0, -a-e)  +  Q2cos  (0+0,+P-e)) 

+(P/2)  (Q,2  +  2  Q,Q2  cos  (20,+a+P)  +  Q/) 

-A  (P/2)  (Q,2cos  (20,+2a+e)  +  2Q,Q2  cos  (a-p+e)  +  Q22  cos  (20,+2p-e)) 


On  substituting 

[A.47] 

0'  -0  +  6,+  a 

[A.  48] 

5  =■  20,  +  2a  +  e 

[A.49] 

<j>  =  a  -  p  +  e 

[A.50] 

the  expression  XX  in  the  joint  probability  density  function  for  r  and  6' 
becomes 
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XX  =  r2  (1-A  cos  (20' -8)) 

-  (2P)as  r  (Q,  COS  8'  +  Q2  cos  (0,-5-h}>)-A  Q,  cos  (0'-8)  -A  Q2  cos  (0'-0)) 

+(P/2)  (0,2(1  -A  cos  5)  +  2  Q,Q2  (cos  (5-<)>)-A  cos  0)  +  Q22  (1-A  cos  (8-20))) 

[A. 51] 

Now  the  signal  phase  0,  is  essentially  random  and  so  to  obtain  the 
probability  density  function  of  r  and  0',  averaged  over  all  possible  signal 
phases,  f(r,0')  is  integrated  with  respect  to  5  to  obtain 

r  exp  ( - ^ - )  I  (x/a^O-A2)) 

f(r,0')  = _ 2og_n-A2]  0 _  [A.52] 

2k  a2  (1  -  A2)0  5 

where  here 

XX=r^  -  (2P)0  Sr  (Q,  cos  0'  -A  Q2cos  (0'-<»)+(P/2)  (Q,2+  2AQ1Q2  cos  <f>  +  Q22) 

[A.53] 

x2=  A2  r4  +  2  A  r3  (2P)0  S  {Q2  cos  (0'-0)-AQ,  cos  0'} 

+r*  P  {Q22(2+A2cos  2(0'-0)+A2Q1 2(2+cos  20')-2AQ1Q2(2cos  0  -cos  (20’-0)} 
-r(2p3)0.5  {A^^cos  0'-AQ12Q2(cos  (8'+0)+2cos  (0'-0))+ 

Q,Q22(2  cos  0'+A2cos  (8’-20))-A  Q23cos  (0'-<|>)} 

+(P2/4)  { A2Q, 4-4 AQ, 3Q2cos  0  +2  Q,2Q22  (2+A2C0S  2p)-4AQ1Q23cos  0  +A2Q24} 

[A.  54] 

If  Q2  =  0  (which  is  approximately  true  in  ail  cases  studied), 

r  exp  (  -if  -.rian0-.5  Q,  CO?  £  t  (E/2)  Q,2)  i  (X  ^(i  -a2)  ) 

f(r,0')  =  _ _ 2  cr2  fl  -  A2  j _ °_  [A.55] 

2k  ct2  (1  -  A2)0  5 

where 

x2  =  A2!4  -  2A2  r3  (2P)0'5  Q,  cos  0'  +  A  r2  P  Q,2  (2+cos  20') 

-  A2  r  (2P3)0  5  Q,3  cos  0'  +  A2  (P/2)2  Q,4  •  [A.56] 

If  A  =  0  (which  is  again  approximately  true  in  the  cases  considered), 
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r  exp  ( -  d:  r  f2P)0,5  a  cos  e1  ±  (P/n  n  2. 
f(r.e')  * _  _ 2ji2_ _ 

2ji  o2 

In  any  case  the  probability  of  detection 

00  2it 

Pd  -  L  /  f(r,0‘)  de’  dr. 

R  0 

In  the  case  where  Q2  =  A  =  0, 

r 

pd  =  r  exPHSr2  +  P  Q,2)/^))  |Q(r  (2Pp  Q1  /(2a2))  dr 
C 00 

e'(u+P)  |q(2V(uP))  du 
J-ln(pfa) 

on  making  the  substitutions 

P  =  PQ,2^) 
u  =  r2/(2o2) 

and  using  equation  [A.33], 


[A.57] 


fA.58] 


[A. 59] 


[A. 60] 


[A.61] 

[A.62] 
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ANNEXE 

ASSUMPTIONS  MADE  IN  CAL  CULATING  PROBABILITIES 

In  deriving  equation  [2.8]  for  the  probability  detection  as  a  function  of 
probability  of  false  alarm  and  the  quantity  P  (defined  by  equation  [2.9]), 
the  following  assumptions  were  made: 

(1)  The  quantities  B  and  C  (see  equations  [2.4]  and  [2.5])  are  normally 
distributed. 

(2)  With  interference  alone,  B  and  C  have  zero  mean. 

(3)  The  quantity  A  (see  equation  [A.1 8])  is  zero. 

(4)  The  quantity  Q2  (see  equations  [A.39]  and  [A.40])  is  zero. 

(5)  The  signal  is  coherent  ie  var(B)  and  var(C)  are  the  same  whether  the 
signal  is  present  or  not. 

Some  of  the  effects  of  relaxing  these  assumptions  will  be  considered 
in  the  remainder  of  this  Annex. 
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B.l  The  Distribution  of  B  and  C 

Suppose  that  B  and  C  are  not  normally  distributed.  Instead  suppose 
that  they  are  independent  (assume  A-0  so  they  have  zero  covariance),  they 
have  zero  mean  and  they  are  identically  distributed  with  probability 
density  function  f.,()  which  depends  only  on  the  variance  cr2  and  not  on  any 

higher  order  moments.  The  probability  density  function  may  then  be 
written 

f,(z)  -f(z/o)/o  [B.l] 

where  z  represents  either  B  or  C  and  the  function  f( )  does  not  depend  on 
any  parameters.  The  probability  of  false  alarm  may  be  derived  as  in  Annex 
A  to  be 


pfa  -J‘  ;r  r  f,  (r  cos  0)  fn(r  sin  0)  d0  dr 
—  2* 

« J  J  u  f(u  cos  0)  f(u  sin  0)  d0  du 
R/o  0 

where  R  is  the  threshold.  Thus  pfa  is  a  monotonic  decreasing  function  of 
R/o  which  has  an  inverse  which  will  be  denoted  by  g(p,a)  -  R/o. 

Similarly,  the  probability  of  detection 

f”  r2* 

Pd  J  rf^rcos  0-B)  f,(r  sin  0-C)  d0dr  [B.4] 

R  0 

••  2k 

»J  s  u  f(u  cos  0-B/o)  f(u  sin  0-C/o)  d0du  [B.5] 

R/a  0 

“  2k 

■  J  /  u  f(u  cos  Q-B  /a)  f(u  sin  0-C  /o)  d0  du  [B.6] 

g(pfa)  0 


(B  2] 
P-3] 


where  B  and  C  are  the  mean  values  of  B  and  C  respectively.  If  f(x)  is  an 
even  function  with  a  single  maximum  at  x«0,  pd  will  be  an  increasing 

function  of  |B  |/o  or  |C  |/o.  For  a  fixed  ratio  |S/C  |,  pd  is  therefore  an 
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* 


) 

I 


< 


f 


increasing  function  of  P  (defined  by  equation  [2.9])  if  Q2  is  approximately 
zero.  Hence  P  is  still  a  sensible  measure  of  detection  performance,  but ,  as 
in  the  case  of  the  normal  distribution  with  A*0,  the  probability  of 
detection  may  depend  on  |S/C  |,  which  in  turn  depends  on  the  signal 
phase.A  measure  of  probability  of  detection  independent  of  phase  may  be 
obtained  (as  in  Annex  A)  by  integrating  equation  [B.6]  with  respect  to 
phase.  Figures  23  and  24  show  the  variation  of  this  integrated  pd  with  P 

for  pfa  between  1 0‘7  and  1 0'2  and 


f(x) 

-0.5exp(-|x|) 

(figure  23) 

[B.7] 

f(x) 

-  (4k)-0  5  { (.8)‘°'S  exp(-1 .25  x2)  + 

(3.2)'0  5  exp(-.3125  x2) } 

(figure  24) 

[B.8] 

The  latter  function  is  the  sum  of  two  normal  distribution  functions. 
Comparison  with  figure  1  shows  that  the  value  of  P  required  to  obtain  a 
particular  probability  of  detection  is  higher  using  [B.7]  or  [B.8]  than  using 
the  gaussian  assumption.  The  difference  for  a  50%  probability  of  detection 
is  5  to  10  dB  using  equation  [B.7]  or  1  to  2  dB  using  equation  [B.8].  The 
difference  is  due  to  the  increased  height  of  the  "tail"  in  equations  [B.7]  and 
[B.8]  (ie  increased  values  of  f(x)  for  large  x)  and  the  consequent  necessity 
of  using  a  larger  threshold  to  achieve  the  same  false  alarm  rate. 
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B.2  The  Mean  of  B  and  C 

If  B  or  C  have  non-zero  mean  (caused  (say)  by  some  dominant 
scatterer  which  is  not  the  target),  then  by  direct  analogy  with  the 
calculations  for  pd  in  Annex  A,  p)#  may  be  written 


Pc- 


e‘(u+Z)  l0(2V(uZ))  du 

R/o 


where  R  is  the  threshhold  and 
Z  m(B2  +  C2)/(2a2) 

with  B  and  C  being  here  the  mean  values  of  B  and  C  with  interference 
alone.  For  a  given  probability  of  false  alarm,  the  threshhold  R  and  hence 
the  probability  of  detection  are  increasing  functions  of  Z.  Figure  25  shows 
the  variation  of  pd  with  P  for  various  values  of  Z  and  pfa  -  10'5.  The  effect 

of  increasing  Z  is  to  increase  the  value  of  P  required  for  a  given 
probability  of  detection.  Even  for  Z-5  (which  would  increase  the  output 
interference  power  by  a  factor  of  6)  the  increase  in  the  value  of  P  required 
tor  a  50%  probability  of  detection  is  only  about  4  dB. 
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B.3  The  Quantity  A 

In  figure  1 .  it  is  assumed  that  A  (defined  by  equation  [A.  18])  is  zero. 
If  A  is  non-zero,  equation  [A.2S]  must  be  used  forf(r)  in  equation  [A. 31]  to 
calculate  the  probability  of  false  alarm.  Further  equation  [A.55]  must  be 
used  for  f(r.e')  in  equation  [A.56]  to  calculate  the  probability  of  detection. 
For  p(>  -  10'5,  figure  26  shows  the  variation  of  probability  of  detection 

with  P  for  various  values  of  A.  From  the  figure,  it  may  be  noted  that 
increasing  A  increases  the  value  of  P  required  for  a  given  probability  of 
detection,  but  this  increase  is  less  than  3  dB  for  the  largest  possible  value 
of  A  (-1 )  for  a  50%  probability  of  detection. 

In  the  cases  considered  in  section  4.1  and  4.2.  where  it  is  assumed 
that  the  carrier  frequency  -  pulse  length  product  is  about  1 5000  (which  is 
fairly  typical).  A  is  always  small  (about  10'2)  and  its  effect  is 
insignificant. 


70 


UNCLASSIFIED 


PROBABILITY  OF  DETECTION 


Figure  2 6 ( U )  Transition  Curves 
for  Selected  Values  of  A 
PF  A  is  1 0  *  * (  -  5 ) 
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B.4  The  Quantity  CL 

In  figure  1 ,  it  is  assumed  that  Q2  (defined  by  equations  [A.39]  and 
[A.40])  is  zero.  If  Q2  is  non-zero,  equations  [A.46]and  [A.51]  may  be  used  in 
equation  [A.58]  to  calculate  the  probability  of  detection.  If  A-0, 

XX  -  r2-(2P)0  5r(Q1cos  0’+Q2cos  (S' -5+4))  +  (P/2) (0,2+20,02008  (5-0)+Q22) 

IB.11] 

in  equation  [A.46],  so  that  on  integrating  equation  [A.58]  with  respect  to  6', 
pd  may  be  expressed  as  in  equation  [A.  59],  but  with  Q,2  replaced  by 

Q,2+2Q,Q2cos  (5-<s)+Q22.  Hence 


Pd“ 


j  e‘(u+p>  l0(2V(uP'))  du 
J-ln(pfa) 


[B.12] 


where 


P'  -  P  (1  +  2kcos  (8-0)  +  k2  )  [B.13] 

k-Q^Q,.  [B.14] 

For  pfa  -  1 0'5,  figure  27  shows  the  variation  of  probability  of  detection 
with  P  for  various  values  of  Q2  averaged  over  the  arbitary  phase  angle 
(5-0)  (with  A-0).  From  the  figure  it  may  be  noted  that  the  effect  of  Q2  is 
minimal  unless  02  is  quite  large  (more  than  0.1).  Increasing  Q2  reduces  the 

value  of  P  required  for  a  50%  probability  of  detection  but  increases  the 
value  required  for  a  high  (eg  90%)  probability  of  detection. 


In  the  cases  considered  in  section  4.1  and  4.2,  where  it  is  assumed 
that  the  carrier  frequency  -  pulse  length  product  is  about  15000  (which  is 
fairly  typical),  CyQ,  is  always  small  (about  10'5)  and  its  effect  is 

insignificant. 
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uiNLohcsoir  ic-u 


Figure  2 7 ( U )  Transition  Curves 
for  Selected  Values  of  Q2/Q1 
PFA  is  1 0*  * ( -5 ) 
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B.5  Signal  Coherence 

Suppose  that,  due  to  multipath  or  other  effects,  the  signal  has  a 
component  which  has  the  same  characteristics  as  interference  ie  the 
contributions  to  B  and  C  from  this  component  have  zero  mean  and  non-zero 
variance.  If  k  is  the  proportion  of  coherent  power  in  the  signal,  then  E[B] 
and  E[C]  are  reduced  by  the  proportion  k.  The  increase  in  the  variance  of  B 
and  C  depends  on  the  exact  nature  of  the  interference-like  component.  If 
this  component  consists  of  a  large  number  of  "signals"  with  random  phase, 
then  it  is  of  the  form 


x(t)  «J*g(e*)  W(t)  cos  (M(t)  +  2icft  +  e1)  de*  [B.15] 

0 

where 

E[g(0')]=O  •  [B.16] 

E[g(e')g(0")]  =  2(i-k)  p  w-eyn.  [B.17] 

The  expected  value  of  x(t)x(u)  is  then 

E[  x(t)x(u)  ]  -  (1-k)  P  W(t)  W(u)  cos  (M(t)-M(u)+2rcf(t-u))  [B.18] 


and  the  additional  variance  in  B  or  C  due  to  the  presence  of  this  extra 
interference  is  (using  equation  [A.1]) 

fV 

s2  -  (1-k)  P  J0  J0  W(t)  W(u)  H(t)  H(u)  cos  (M(t)-M(u)+27tf(t-u)). 

cos  (N(t)-N(u)+2rrf(t-u))  dt  du  [B.  1 9] 

(V 

-  0.5  (1-k)  P  JQ  JQ  W(t)  W(u)  H(t)  H(u)  cos  (M(t)-M(u)-N(t)+N(u))  dt  du 

[B.20] 

-  (1-k)  P  0,2/2  [B.21] 

-  2  (1-k)  Po2  [B.22] 

Now  the  probability  of  false  alarm  may  be  calculated  using  equation  [A. 33] 
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in  its  present  form,  but  the  probability  of  detection  is  calculated  from 
equation  [A.59]  with  o2  replaced  by  (  o2+  s2).  Hence 


Pd  = 


e-(u+z)  |q(2V(uZ))  du 

-a2  ln(pf gj/fa2  +S^) 


where 

P(k)  =  kPQ^/^As2) 

=  k  P(0)  /  (1  +  2(1-k)P(0) ) 

and  P(0)  =  P  C^2  /  4 o2  is  the  value  of  P  (defined  by  equation  [2.9])  if  k  were 
zero.  For  pfa  =»  1 0'5,  figure  28  shows  the  variation  of  probability  of 

detection  as  a  function  of  P(0)  for  various  values  of  k.  It  may  be  noted 
from  the  figure  that  decreasing  k  will  decrease  (slightly)  the  value  of  P 
required  for  a  50%  probability  of  detection,  but  will  increase  (by  up  to  6 
dB)  the  value  of  P  required  for  a  90%  probability  of  detection. 
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igure  2 8 ( U )  Transition  Cunves 
fon  Selected  Values  of  k 
PFA  is  1 0  *  * ( -  5 ) 
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ANNEX C 

ASSUMPTIONS  MADE  IN  CALCULATING  GAINS 

In  this  Annex,  generalised  versions  of  equations  [4.12]  and  [4.14]  will 
be  derived  (from  which  these  equations  follow  as  a  special  case).  These 
equations  and  the  results  of  Annex  A  will  be  used  to  examine  the  effect  of 
modifying  some  of  the  assumptions  made  in  sections  4  and  5.  These 
assumptions  are  listed  at  the  end  of  this  part  and  the  remaining  parts 
discuss  the  modification  of  each  assumption  in  turn. 

Suppose  that  the  covariance  of  g  (see  section  4  for  a  definition  of  g) 
is  of  the  form: 

E[g(t,C\f,e')  g(u,?",f",6")  ]  -  2  n  V(t-u)  L(f)  5(f-f")  X(C-0  8(6'-9")/(nK) 

[C.1] 


where  the  amplitudes  g  are  independent  for  different  phases  and  centre 
frequencies,  L(f)  gives  the  frequency  distribution  of  the  reverberation 
power,  X  gives  the  correlation  with  different  path  lengths  (a  delta 
function  in  section  4),  V  gives  the  correlation  at  different  times  (constant 
in  section  4)  and  the  functions  are  scaled  so  that  the  average 
reverberation  power  over  the  the  analysis  period  is  R.  The  functions  V  and 
X  are  even.  The  expected  value  of  x(t)x(u)  may  be  written 

E[  x(t)x(u)  ]-  (R/K)  V(t-u)  J  *  /  **  L(f')  X(C-0  W(K')  W(u-?"). 

t-K  u-K 

f  cos(M(t-Q+2nft+6')  cos  (M(u-?’,)+27rfu+e')  d0'  df  d^"  d£'  [C.2] 

=  (R/K)  V(t-u)  /*  JU  1°°  L(f)  X(C’-C)  W(t :<)  W(u-0- 

t-K  u-K  -o° 

cos  (M(K')-M(u-0+2rcf(t-u))  df  dC'  dC  [C.3] 
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-  (R/K)  V(t-u)  { L(f)  cos  27tf (t-u)  d f  . 
f'  ru 

^  W(t<r>  W(u'W  cos  (M(t-C)-M(u-O)  dC"  df 

-  L(f)  sin  27tf(t-u)  df  . 
ru 

W<u'^")  sin  (M(t-C)-M(u-O)  dC  df 
Now,  on  letting  x  =  u-f  and  y  -  t-f  and  v=t-u, 
rf  ru 

|.KJU-KX(^')  W(t'^  W^-^)  cos  (M(t-f)-M(u-f))  dr  df 

mOo  X(v+x*y)  WM  cos  (M(y)-M(x))  dy  dx 
*  W(v)  (say). 

W(v)  is  an  even  function  of  v  (as  may  be  seen  by  interchanging  the  order  of 

ml  h  h  n0ti"9  that  X  iS  eVen)'  An  a,tematlve  expression  for  W(v) 
Tin3  ^ in"erinteg'a(  into  Parts:  from  o  to 

coZrng71ZC^0Z9^  °rd6r  °f  int6grati0n  in  ^ firSt  P3rt  and 

W(v)  =  q  /x  (X(v+x-y)+X(v-x+y))  W(y)  W(x)  cos  (M(y)-M(x))  dy  dx  i 

K  K-x  ^ 

3  {  {  (x(v+u)+X(v-u))  W(u+x)  W(x)  cos  (M(u+x)-M(x))  du  dx  [ 

"  o  Jo  (x(v+u)+x(v-u))  W(u+x)  W(x)  cos  (M(u+x)-M(x))  dx  du  [< 

on  letting  u  =  x-y  and  changing  the  order  of  integration  On  makino  a 
further  substitution  w  =  K-2x-v  in  the  inner  integral 


W(v)  _  K+u(X(v+u)+X(v-u))  W,(x-u)  W,(u+x) 


where 


wi(x)  =  W(  (K-x)/2 ) 


cos(M,  (x-u)-M,  (u+x))  dx  du 

[C.10] 

[C.11] 
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M,(x)  -M((K-x)/2).  [C.12] 

Expressions  for  the  second  double  integral  in  equation  [C.4]  in  which  sin 
replaces  cos  may  be  derived  in  a  similar  way  to  obtain 

W, (v)  -  /  JU  X(C'-C)  W(t-C)  W(u-?")  sin  (M(t-C)-M(u-C"))  d£"  d£  [C.13] 

t-K  u-K 

-  JK/K"U(X(v+u)-X(v-u))  W(u+x)  W(x)  sin  (M(u+x)-M(x))  dx  du  [C.14] 

0  0 


0.5]*}  (X(v+u)-X(v-u))  W,(x-u)  W,(u+x)  sin(M1(x-u)-M1(u+x))dxdu 


where ' 


0  -K+u 

is  an  odd  function. 


[C.15] 


If  W  and  M  are  symmetric  about  K/2  (which  is  true  for  the  CW  and  FM 
chirp  pulses  and  the  common  shading  functions),  W,  and  M1  are  even 

functions  and 


K  K-u 

W(v)  *  j  }  (X(v+u)+X(v-u))  W,(u-x)  W^u+x)  cos  (M^u-xJ-M^u+x))  dx  du 

[C.16] 

W^v)  m  0.  [C.17] 


Additionally,  if  X  is  the  Kroneker  delta  function,  as  in  sections  4  and  5, 

K-v 

W(v)=J  W(v+x)  W(x)  cos  (M(v+x)-M(x))  dx  [C.18] 

o 

K-v 

-  J  W^v-x)  W^v+x)  cos  (M^v-xJ-M^v+x))  dx.  [C.19] 

Now  consider  the  other  integrals  in  equation  [C.4],  If  L(f)  has  a  main 
peak  at  f«p,  define 

on 

L(v)  -  V(v)  |  L(x+p)  cos  2rcfxv  dx  [C.20] 

-oo 
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Mv)  -V(v)  /  L(x+p)  Sin  2tzxv  dx. 

-•o 

The  integrals  in  equation  [C.4]  may  then  be  written 

V(t-U)  L(f)  cos  2nf(t-u)  dr.  L(t-u)cos  2w(t-u)  -  L,(t-d)sin  a^fru) 

VM  CL<f> si"  «f-  MMsin  avft-u)  ♦  L,(t-u)cos  2*Kt-u] 


If  L (f)  is  symmetric  about  f=|i,  as  occurs  with  the  functions 
sections  4  and  5, 


considered  in 


fC  2i] 


[C.22] 


[C.23] 


L(v)  =  2  V(v)  L(x+|a)  cos  2nxv  dx 
L,(v)  =0. 

In  any  case, 


[C.24] 


[C  .25] 


E[  x(t)x(u)  I  -  R  { ( L(t-u)  W(t-u)  -  L,(t-u)  Wt(t-u) )  cos  2*p(t-u) 

-  (  L|  (t-u)  W(t-u)  +  L(t-u)  W,(t-u) )  sin  2nyi(i-u)  }.[C.26] 

which  is  an  even  function  of  (t-u)  only  and  will  be  denoted  Z(t-u)  The 
average  reverberation  power  over  the  analysis  period  is  R  if 

1(0)  W(0)  =  K. 

bV  ™  50  «  eQuation  [C.27]  scales 

In  the  cases  considered  in  section  4  and  5,  equation  [C.26]  becomes 
E[  x(t)x(u)  ]  *  (R/K)  L(t-u)  W(t-u)  cos  27tp(t-u) 

where  L  and  W  are  given  by  equations  [C.24]  and  [C.18]  with  V(v)=1  From 
equations  (2.2]  and  (C.iaj.  W(0)-K  so  that  equation  (C.27,  imp, LI 


[C.27] 


[C.28] 
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L(0)-1.  [C.29] 

Substituting  Z(t-u)  for  E[  x(t)x(u)  ]  in  equation  [2.1 1]  gives 
fT  fT 

o2  -  J0  Z(v)  Jv  H(t)  H(t-v)  cos  (N(t)-N(t-v)+2itfv)  dt  dv  [C.30] 

J  rT'v 

"0  5  o  j  H,(w-v)  H,(w+v)  cos  (N1(w-v)-N1(w+v)+2jcfv)  dt  dv 

[C31] 

where 


^(x)  -H((T-x)/2)  [C.32] 

N,(x)  -  N{  (T-x)/2  ).  [C.33] 

Define 

H(v)  -  J  H(t)  H(t-v)  COS  (N(t)-N(t-v))  dt  [C.34] 

V 

.T-v 

-  0.5  J  ^  H,  (w-v)  H,  (w+v)  cos  (N1  (w-v)-N,  (w+v))  dw  [C.35] 

T 

H^v)  -  J  H(t)  H(t-v)  sin  (N(t)-N(t-v))  dt  [C.36] 

T-v 

-  0.5  J  H^w-v)  H^w+v)  sin  (N,  (w-v)-N,  (w+v))  dt  [C.37] 


which  are  directly  analogous  to  the  definitions  of  W  and  W, .  Equation 
[C.31]  then  becomes 

J 

o2  -  (R/K)  J0  { (  L(v)  W(v)  -  L,  (v)  W,  (v ) )  H(v)  cos  2rcpv  cos  2itfv 

-  (  L^v)  W(v)  +  L(v)  W,(v  ) )  H(v)  sin  2k\l v  cos  2?tfv 

-  (  L(v)  W(v)  -  L,(v)  W,  {v  ) )  H,(v)  cos  2k\iv  sin  2rcfv 

+  (  L,(v)  W(v)  +  L(v)  W,(v  ) )  H,(v)  sin  27t|iv  sin  2irfv  }  dv. 

[C.38] 

Each  of  the  products  of  sin  or  cos  may  be  written  as  the  sum  of  two  sin  or 
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cos  terms  with  argument  2jt(f-p)v  and  2ji(f+p)v.  The  integrals  containing 
the  second  of  these  arguments  will  be  small  compared  with  those 
containing  the  first  if  (f-p)T  «  (f+p.)T,  as  in  practical  situations.  Thus  in 
the  cases  considered  in  section  4  and  5  where  H  and  N  are  symmetric  about 
T/2, 

T 

a2  -  •q  L(v)  W(v)  H(v)  cos  2jr(f-p)v  dv  [C.39] 


H(v)  -  H,(w-v)  H,(w+v)  cos  (N,(w-v)-Nt(w+v))  dw  [C.40] 

H,(v)  -0.  [C.41] 


The  additional  assumptions  made  in  sections  4  and  5  are  listed  below. 
Some  effects  of  varying  these  assumptions  are  described  in  the  following 
sections. 

(1 )  The  start  of  the  signal  coincides  with  the  start  of  the  analysis  period 
ie  D=0  in  equation  [2.7], 

(2)  The  analysis  and  signal  frequencies  coincide  ie  f,*f. 

(3)  The  duration  of  the  signal  equals  the  duration  of  the  analysis  period 


(4)  The  scatterers  have  constant  properties  for  times  equal  to  the  length 
of  the  analysis  period  ie  V(v)»1 . 

(5)  The  scatterers  are  independent  for  different  path  lengths  ie  X(v)»8(v). 
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C.1  The  Start  of  the  Signal  and  Analysis  Period  Coincide 

The  condition  D-0  (ie  the  start  of  the  signal  and  analysis  periods 
coincide)  is  required  for  equations  14.20]  or  [4.27]  to  be  valid  for  Q1 . 

However  if  D>0  (but  all  the  other  assumptions  are  unchanged),  the  values 
of  o2  are  unchanged  but  is  given  by 

rT 

Q2  »  (JQ  W(t-D)  H(t)  COS  (  M(t-O)  -  N(t)  )dt)2  + 

(J0  W(t-D)  H(t)  sin  (  M(t-D)  -  N(t)  )dt)2  [C.42] 

If  D<0  a  similar  expression  for  Q,2  may  be  written  with  different 
integration  limits.  The  lower  limit  becomes  zero  and  the  upper  limit 
T+D-T- 1 D I .  On  letting  t*«t-D,  an  equation  equivalent  to  [C.42]  may  be 
derived  in  which  H  and  W  are  interchanged,  M  and  N  are  interchanged  and  D 
is  replaced  by  j  D  j .  In  the  following  D>0  only  will  be  considered,  as  D<0 
gives  identical  results  by  interchanging  shaping  and  shading. 

For  CW  pulses  and  each  of  the  combinations  of  H(t)  and  W(t)  used  in 
section  4,  the  degradation  in  P  is  presented  in  figure  29  for  values  of  D/T 
between  0  and  1 .  This  is  the  same  as  figure  2  except  that  here  the 
degradation  is  taken  relative  to  its  value  when  D-0,  not  necessarily 
relative  to  its  optimal  vaiue.lf  there  is  no  overlapping  of  analysis  periods, 
the  maximum  degradation  will  occur  if  the  signal  starts  at  time  T/2  (if  it 
starts  at  a  later  time  during  the  period,  the  next  period  will  contain  the 
signal  for  a  longer  time  and  it  should  be  considered  with  D<0,  giving 
identical  results  to  those  in  the  figure)  and  corresponds  to  D/T -0.5.  If  the 
analysis  periods  overlap,  the  maximum  will  correspond  to  a  lower  value  of 
D/T,  but  more  processing  will  be  required  and  the  false  alarm  rate  may  be 
increased  due  to  the  increase  in  the  number  of  analysis  periods  per  unit  of 
time.  Of  course  the  overlapping  periods  will  not  be  independent  so  that 
this  increase  may  not  be  significant. 

For  FM  pulses  with  M(t)-N(t), 
fT 

Q,2  -  (JQ  W(t-D)  H(t)  cos  (  M(t-D)  -  M(t)  )dt)2  + 

J 

(}D  W(t-D)  H(t)  sin  (  M(t-D)  -  M(t)  )dt)2  [C.43] 
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For  a  chirp  pulse  where  M(t)  is  given  by  equation  [4.23], 
fT 

(V  -  (JD  W(t-D)  H(t)  cos  2np  {  2D(t-T/2)  -  D2  }  dt)2  + 
rT 

(JD  W(t-D)  H(t)  sin  2rrp  {  2D(t-T/2)  -  D2 }  dt)2. 

-  T2(^d/t  W((u-D/T)T)  H(uT)  cos  (  2rcDu/S  )du)2  + 

T2(JD/t  W((u-D/T)T)  H(uT)  sin  (  2k  Du/S  )du)2  [C.44] 

where  S-(2pT)'1  is  the  inverse  of  the  FM  bandwidth.  When  using  FM  pulses, 
the  analysis  periods  overlap  and  the  start  of  successive  intervals  differ  by 
a  time  AT  which  should  be  set  less  than  or  equal  to  S.  Figure  30  shows  the 
degradation  in  the  gain  for  values  of  D/S  between  0  and  5  (assuming  that 
S/T  «  1  ie  the  bandwidth-  pulse  length  product  is  large  so  that  the  lower 
limits  of  integration  in  equation  [C.44]  may  be  set  equal  to  zero).  If  the 
starting  times  of  the  analysis  periods  are  separated  by  aT-S,  the 
maximum  degradation  will  occur  when  D/S  -  0.5  and  is  about  4  dB  with  no 
shading  or  shaping,  but  is  rather  less  (at  most  1 .5  dB)  otherwise.  For  small 
values  of  D/S  (less  than  1),  the  degradation  is  greatest  when  there  is  no 
shading  or  shaping,  indicating  that  a  plot  of  processor  output  as  a  function 
of  analysis  period  start  time  would  have  a  sharper  peak  corresponding  to 
the  signal  in  this  case.  However,  for  larger  values  of  D/S  (at  least  3)  the 
situation  is  reversed  and  the  degradation  is  greater  when  the  shading  or 
shaping  schemes  are  used.  Thus  the  "side  lobes"  in  the  plot  will  be  lower. 

This  lowering  of  side  lobes  has  important  implications  if  there  is  another 
unwanted  signal.  The  use  of  shading  and  shaping  may  reduce  its 
interference  effect,  even  if  it  is  quite  large. 
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Figure  2S(U)  Degradation  in  Gain 
as  Funct i on  of  D/1  for  CW  Pulse 
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DEGRADATION  IN  P  (DB) 
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C.2  The  Analysis  and  Signal  Frequencies  Coincide 

The  condition  f-f,  (ie  the  centre  analysis  frequency  is  the  same  as 

the  centre  signal  frequency)  is  required  for  equations  [4.20]  or  [4.27]  to  be 
valid.  However,  if  f^f,  (and  all  other  assumptions  are  unchanged  except 

that  possibly  D>0),  the  values  of  a2  are  unchanged  but  Q,  is  given  by 
J 

Q,2  -  (JD  W(t-D)  H(t)  cos  (  M(t-D)  -  N(t)  +  2n(f-  f,)t  )dt)2  + 
fT 

(JD  W(t-D)  H(t)  sin  (  M(t-D)  -  N(t)  +  2n(f-  f,  )t  )dt)2  [C 

For  a  CW  pulse, 
fi 

Qi2  -  t2(Jd n  W((u-D/T)T)  H(uT)  cos  (  2*  If-  f,  I  Tu  )du)2  + 

T2(Jd/t  W((u-D/T)T)  H(uT)  sin  (  2rc  |  f-  f,  I  Tu  )du)2  [C 

which  is  exactly  the  same  as  equation  [C.44]  with  D/S  replaced  by  (f-f^T, 
so  that  figure  30  also  shows  the  degradation  in  P  due  to  the  frequency 
mismatch  for  values  of  I  f-fn  |  T  between  0  and  5  if  D=0.  As  it  would  be 

expected  that  the  frequency  bins  would  be  separated  by  at  most  Af=T '1 , 
the  maximum  loss  due  to  frequency  mismatch  would  correspond  to 
I  f-f1 1  T-0.5  or  4  dB  with  no  shading  or  shaping.  As  for  the  FM  case,  the 
signal  spectral  peak  is  broader  using  shading  or  shaping,  but  the  side  lobes 
are  lower,  improving  the  rejection  of  a  strong  unwanted  signal  at  another 
frequency. 

The  combined  effect  of  having  D*0  and  tef,  is  illustrated  in  figure 
31 ,  where  the  degradation  in  P  is  plotted  as  a  function  of  (f-f,)T  for 
selected  values  of  D/T  and  no  shading  or  shaping.  Notice  that  the  worst 
case  of  D/T-0.5  and  (f-f,)T»0.5  gives  a  degradation  of  about  6  dB. 

For  an  FM  chirp  pulse, 


87 


UNCLASSIFIED 


UNCLASSIFIED 


Qi2  *  t2(JD/t  W((u-D/T)T)  H(uT)  cos  (  2re(D/S+(f-  f^TJu  )du)2  + 

T2(jD/t  W((u-D/T)T)  H{uT)  sin  (  2n(D/S+(f-  f^JJu  )du)2 

In  this  case,  the  detection  algorithm  is  performed  at  one  centre  frequency 
only  (so  that  f-F,  the  transmitted  pulse  centre  frequency).  Hence,  even  if 
D-0,  I  f-f,  I T  can  be  quite  large  if  there  is  significant  doppler  shift  in  the 

signal  return.  If  I  f-f  1 1  T>0.5,  a  different  analysis  period  would  give  the 
largest  processor  output  so  the  estimated  range  to  the  target  would  be  in 
error  -  the  largest  output  may  be  estimated  by  adding  an  integer  multiple 
of  S  to  D  such  that  q*  I  D/S+(f-f1)T  I  <0.5.  Figure  30  shows  the  degradation 

in  P  as  a  function  of  q  in  the  case  where  1 D  |  /T  is  zero  or  close  enough  to 
zero  to  be  ignored.  The  degradation  for  values  of  q>0.5  correspond  to  the 
reduction  in  echo  in  range  cells  neighbouring  the  one  with  maximum 
response.  If  the  frequency  shift  is  such  that  |D  j/T  may  not  be  ignored, 
figure  31  shows  the  degradation  in  P  as  a  function  of  q  and  selected  values 
of  |  D  |/T  in  the  case  of  no  shading  or  shaping. 


[C  47] 
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C.3  The  Pulse  and  the  Analysis  Period  Have  Equal  Length 

Equality  in  the  length  of  the  pulse  K  and  the  length  of  the  analysis 
period  T  is  necessary  for  the  equation  [4.20]  or  [4.27]  for  Q.,2  to  be  valid. 

If  T*K,  these  equations  should  be  replaced  by 

jnin(T,(T+KV2] 

Qi2  -  (W(T-KV2]  W(t-(T-K)/2)  H(t)  dt)2  [C.48] 

if  it  is  assumed  that  the  signal  pulse  and  analysis  period  are  aligned  so 
that  their  centres  coincide  (Ie  the  frequency  modulation  of  the  signal  and 
the  analysis  replica  are  the  same  throughout  the  integration  period).  This 
alignment  gives  the  gretest  value  of  Q,2.  The  values  of  cr2  are  also  changed 

in  this  case.  Since  W(v)-0  unless  0£v£K, 

rinfT.K] 

0  L(v)  W(v)  H(v)  cos  2n(f-|i)v  dv  [C.49] 

v 

-  W/2  j  0  L(Ku)  W(Ku)  H(Ku)  cos  cou  du  [C.50] 

where  u  -  v/K  ,  -  min[1  ,T/K]  and  co  -  2n(f-p)K. 

For  CW  pulses,  the  gain  as  a  function  of  to  is  shown  in  figures  32  to  40 
in  the  case  where  L(v)«1  (case  1  from  table  1)  for  T/K  taking  the  values 
0.1 ,  0.5,  0.9, 1 .0. 1 .1  and  2.0.  The  case  T/K  ■  1 0.  is  also  included  in  figures 
36,  38  and  40,  but  not  in  the  others  as  the  line  is  strongly  oscillatory  and 
tends  to  obscure  the  other  lines.  The  figures  show  the  9  possible 
combinations  of  shading  and  shaping  functions  from  tables  2  and  3.  In  this 
case,  the  quantity  J  from  section  5.1  (the  integral  in  equation  [C.50])  may 
be  written 


J  -  Z(-1)n 

n-0 

where 


f  S^^v)  sin  coy  s^2n+1'(\|/)  cos  coxy  -  S^2n+1*(0) 
-  + - 


V 


co 


,2n+1 


CO 


2n+2 


[C.51] 


S(u)  -W(Ku)H(Ku). 
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Now  S(\y)-0  since  W(K)-0  and  H(T)»0,  but  in  this  case, 


S*(V) 


KW(T)H'(T)  if  T  <  K 
K  W(K)  H{K)  if  T  >  K. 


[0.52] 


For  the  shading  and  shaping  functions  from  tables  2  and  3,  S’(y)  is 
non-zero  forT<K  unless  cos  shading  is  used  and  S'(y)  is  non-zero  for  T>K 
unless  cos  shaping  is  used.  This  is  different  from  the  results  of  section 
5.1  and  explains  the  oscillations  in  figures  32  and  33  and  the  lack  of 
oscillations  for  some  lines  in  figures  37,  38  and  40.  Note  also  that 

S'(O)  -  K  ( T  W(0)  +  K  H'(0) )  [C.53] 

so  that  S'(0)  is  non-zero  unless  both  cos  shading  and  shaping  are  used  as  in 
section  5.1 .  Hence  the  gain  will  increase  with  to2  (ie  6  dB  per  doubling  of  co) 
except  where  cos  shaping  and  shading  are  used.  If  both  cos  shaping  and 
shading  are  used,  S'(0)-S'(v)“S"(v)“Sm(0)"S"'(v)“S<4,(v)“0  and 


S(5)(V) 


S<5>(0) 


-KsW(T)H"(0)  if  T<  K 

.  -K5  W"( 0)  H(K)  if  T  >  K. 

K5(TW<5>(0)  +  KH^(0)) 
-K5(T{W"(0)}2  +  K{H"(0)}2) 


[C.54] 

[C.55] 


Hence  the  gain  oscillates  in  this  case  also  but  is  proportional  to  co6  as  may 
be  seen  in  figure  34.  In  figures  35  and  36,  it  would  be  expected  that  there 
would  be  oscillatory  behaviour  for  both  T<K  and  T>K,  since  S'(y)*0. 
However  S*(y)  is  very  small  if  T<K  in  figure  35  or  if  T>K  in  figure  36  since 
H(0)  or  W(0)  are  very  small  in  the  respective  cases.  Thus  the  oscillations 
are  present  but  not  noticeable  in  these  figures.  In  figure  39,  it  would  be 
expected  that  there  would  be  no  oscillation  if  T<K  since  S'(v)-0.  However, 
S'(O)  is  small  and  for  the  values  of  co  being  considered,  an  oscillation  from 
a  higher  order  term  in  equation  [C.51]  gives  the  (reducing)  oscillatory 
behaviour  in  this  case.  This  effect  also  explains  the  small  oscillations 
then  T/K-0.1  in  figure  37  and  when  T/K-10  in  figure  40. 
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For  an  FM  chirp  pulse,  the  analysis  of  section  5.2  gives  the  result 
(similar  to  equation  [5.25])  that  when  co-0 
fw 

J  -  TV1  J0  L(vK/y)  W(vK/y)  H(vK/y)  dv 


where 


y  -  2jtpK2 
y  =  min  [1  ,T/K]. 


Hence 


J  '  X  T2/y 


where 


X  =»  lim  ^  L(vK/y)  W(vK /y)  H(vK/y)  dv. 


Jgh^v)  h2(v)  dv. 


[C.56] 


[C.57] 

[0.58] 


[C.59] 

[C.60] 

[C.61] 


where 


h,(v)  -  W,2(Kw)  cos  vwdw  [C.62] 

r1 

h2(v)  =  J0  H12(Kw)  cos  vw  dw  [C.63] 

Substituting  equations  [C.62]  and  [C.63]  into  equation  [C.61]  gives 


X  -  ton  %(X) 


where 


X(X)  -  0.5  J0  J0  H12(Kw)  W^fKx)  ^0(cos(w-x)v+cos(w+x)v)  dvdx  dw 


0.5  ^  H,2(Kw)  W12(Kx)  sin  (w-x)X/  (w-x)  dx  dw  +  QQ 


[C.64] 

[C.65] 


101 


UNCLASSIFIED 


UNCLASSIFIED 


where  the  additional  terms  are  grouped  together  into  the  term  QQ.  These 
terms  QQ  may  all  be  expressed  in  the  form 

.b 

J  g(x)  sin  Xx  dx 

a 

where  g(x)  is  integrate  in  the  Riemann  sense.  Riemann's  lemma(see,  for 
example  reference  [4],  page  40)  states  that  integrals  of  this  form  approach 
zero  as  X-*»,  so  that  QQ->0  as  X->«  and  QQ  .  Ignoring  QQ  and  letting 
u-(w-x)  in  equation  [C.65], 

fv 

*(X)  .  J0  f(u)  sin  Xu  /  u  du  [C.66] 

where 

f(u)  -  Ju  H12(Kw)  W12(K(w-u/X))  dw.  [C.67] 

Provided  that  (f(u)-f(0))/u  is  integrable  for  0<u<\y  (which  will  be  true  for 
any  sensible  shaping  or  shading  functions),  Riemann's  lemma  implies 

J0  (f(u)  -  f(0))  sin  Xu  /  u  du  -»0  [C.68] 

as  X— Hence 


X  -  f(0)  lim  J0  sin  u  /  u  du 

X—*oo 

-  (jt/2)  JQ  H,2(Kw)  W12(Kw)  dw. 

This  justifies  equation  [5.33]  in  the  case  where  T-K. 

By  substituting  the  expression  for  J  given  by  equations  [C.59]  and 
[C.70]  into  equation  [5.1]  and  noting  that  equation  [C.48]  may  be  rewritten 

rv 

Q,2  -  K2  (J„  W,(Ku)  H,(Ku)  du)2  [C.71] 
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the  gain  G  may  be  written 

G  -  y  (J,  W,(Ku)  H,(Ku)  du)2  /  n  ^  W,2(Ku)  H,2(Ku)  du. 

This  may  be  maximised  by  choosing  shading  and  shaping  functions  such 
that  W^Ku)  H^Ku)  -  1  eg  no  shading  or  shaping.  The  gain  is  independent  of 

T/K  for  T/KS1  if  there  is  no  shading,  but  if  other  shading  functions  are 
used,  since  H,  depends  on  T,  the  gain  does  depend  on  T/K. 
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C.4  The  Scattered  have  Constant  Properties 

This  condition  is  necessary  in  order  that  V(v)«1  and  the  functions  L(v) 
in  table  1  be  valid.  If  the  scatterers  do  not  have  constant  properties  (ie 
the  target  strength  or  speed  or  position  of  the  scatterers  changes  during  a 
time  period  equal  in  length  to  the  pulse),  then  V(v)  will  be  a 
non-increasing,  non-negative  function,  scaled  so  that  V(0)-1.  The 
functions  L(v)  in  table  1  should  then  be  multiplied  by  V(v).  If  the  new 
function  L(v)  is  of  the  form  of  one  of  the  functions  in  table  1 ,  then  the 
graphs  in  figures  4  to  21  may  be  used  to  deteremine  the  gain  for  a  CW 
pulse.  Alternatively,  the  analysis  in  section  5  shows  that  for  large  to  the 
gain  depends  only  on  L'(0)  (or  if  this  is  zero,  on  the  lowest  derivative  of  L 
which  is  non-zero  at  the  origin).  The  gain  for  an  FM  pulse  is  probably 
unaffected  as  previously  noted  in  sections  4  and  5  for  ail  but  small  values 
of  y. 
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C.5  The  Scattered  are  Independent 

If  the  scatterers  are  independent  for  different  path  lengths,  then 
X(v)-S(v)  in  equation  [C.1].  However,  multiple  paths  to  and  from  a  single 
scatterer  may  have  different  lengths,  and  it  is  possible  that  there  may  be 
some  small  correlation  between  scattered  pulses  arriving  at  slightly 
different  times,  so  that  X(v)  may  be  non-zero  over  a  region  of  non-zero 
length.  Alternatively,  the  scatterers  may  form  a  periodic  structure  (eg 
ocean  waves)  so  that  the  amplitudes  of  scattered  pulses  arriving  at  times 
differing  by  a  constant  may  be  correlated.  In  this  case,  X(v)  may  be  of  the 
form 


X(v) 


N 

%  8(v)  +  I  a,  ( 5(v-Vj)  +  Sfv+Vj) ) 


i-i 


Although  the  form  of  X(v)  may  not  be  known,  the  following  gives  an 
indication  of  the  changes  required  to  the  analysis  of  section  5  if  X(v)*5(v). 

The  function  W  (defined  by  equations  [C.5]  and  [C.6])  may  be  written 

W(v)  -  f*  ( *(v+u)  +  X(V'U) )  W0(u)  du 

where  W0(u)  is  the  function  W  in  the  case  where  X(v)  is  equal  to  5(v).  The 
function  X  is  scaled  so  that  W(0)  -  K  as  before.  The  expressions  for  Q^,  1, 
H  that  are  used  in  sections  4  and  5  remain  unchanged. 

For  a  CW  pulse,  the  important  difference  of  a  more  general  X(v)  on  the 
analysis  of  section  5.1  is  that,  even  for  cos  shaping,  W(t)  may  not  be  zero. 

In  this  case,  the  gain  will  be  oscillatory  for  large  w,  since  S,2n+1)(1)  will 
be  non-zero  for  the  smallest  value  of  n  for  which  S<2n+1)(0)  is  non-zero. 

For  example,  if  L(v)-1  and  H’{0)*0,  then 

S'(0)  -  T2{  H’(0)+W(0) )  *  0 

S(1)  -  TH’(T)W(T)  *  0 

since  H'(T)  -  H'(O).  Thus, 
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J  „  T  (  H*(T)  W(T)  cos  co  -  T  (  H*(0)+W(0)) )  /  co2  [C.73] 

and  the  gain  is  proportional  to  co2  and  oscillatory.  If  X(v)  may  be  written 

X(v)  -  a  5(v)  +  b  X,  (v)  [C.74] 

where  X,(v)  is  even,  non-negative  and  is  either  bounded  or  consists  of  the 
sum  of  delta  functions,  none  of  which  is  centred  at  v-0  or  v«±K,  then  the 
contribution  to  W(0)  from  X,(v)  is  zero  and 

W(0)  -  a  W0'(0).  [C.75] 

However, 

W(T)  -  JT  (  X,  (T+u)  +  X,  (T-u) )  W0(u)  du  [C.76] 

so  that  an  increase  in  b,  which  probably  gives  a  decrease  in  a  in  order  that 

W(0)  -  a  K  +  2  b  T  X,(u)  W0(u)  du  -  T,  [C.77] 

increases  the  oscillatory  component  of  the  gain  but  may  decrease  the 
non-oscillatory  component.  If  W(T)  -  0  (which  would  be  fortuitous  rather 
than  planned  if  b*0  due  to  the  uncertain  nature  of  the  function  X(v)),  the 
oscillatory  term  is  zero  and  the  larger  b  is  compared  with  a,  the  larger 
will  be  the  gain.  In  fact  if  a-0  (as  assumed  incidently  in  the  base  case  in 
section  4),  W(0)-0  even  with  a  shaping  function  which  is  non-zero  at  the 
endpoints  of  the  pulse.  In  this  case,  selecting  the  shading  function  such 
that  H(0)-H(T)-0  so  that  H’(0)-0  is  sufficient  to  ensure  that  the  gain  is 
proportional  to  co4  ( or  higher  powers),  still  assuming  that  l(v)-1 .  The 
shaping  function  is  then  only  relevant  in  determining  the  proportionality 
factor. 
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